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HSU\-100  is  a  100ksi  nominal  yield  strength  steel  being  developed  by  the  Navy  for 
Naval  shipbuilding  applications.  To  assist  in  ductile  and  brittle  failure  modelling  of  this  low 
carbon  steel,  tensile  tests  were  conducted  at  temperatures  ranging  from  37°C  to  liquid 
nitrogen  (-196°C)  to  determine  a  constituitive  equation  for  this  alloy  at  low  temperatures. 

The  Hollomon  Power  Equation  (a=  Ken)  and  the  Voce  Equation 

(a=aoo-[aoo-a0]exp[-€/A])  are  used  to  describe  the  true  stress  /  true  strain  behavior 

to  failure  of  individual  tests.  The  Bridgman  correction  is  applied  to  the  true  stress  to 
compensate  for  the  triaxial  stress  state  that  exists  at  and  beyond  the  necking  (maximum 
load)  point.  The  Power  and  Voce  equations  are  then  fit  to  the  Bridgman  corrected  true  stress 
versus  true  plastic  strain.  Relative  comparisons  are  made  between  the  two  equations 
resulting  in  the  conclusion  that  the  Voce  equation  describes  the  stress-strain 
characteristics  of  this  alloy  better  than  the  Hollomon  equation. 

The  temperature  dependence  of  the  material  constants,  oQ,  o^,  and  A,  in  the  Voce 

equation  were  determined  producing  a  constituitive  equation  for  the  tensile  behavior  of 
HSLA-100  as  a  function  of  strain  and  temperature,  at  a  strain  rate  of  about  1x10"3 
in/in-sec.  The  Bridgman  correction  factor  was  investigated  in  detail,  mapping  its  change  as 
a  function  of  strain  past  the  necking  point.  The  Bridgman  correction  produced  a 
discontinuity  in  the  stress-strain  curve  at  strains  in  the  vicinity  of  the  necking  strain. 
This  is  not  believed  to  be  the  true  material  behavior  and  thus,  indicates  that  the  Bridgman 
correction  methodology  is  suspect. 
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I.  INTRODUCTION 

A.   CLASSIFICATION  OF  HSLA  STEELS 

The  term  HSLA,  High-Strength  Low-Alloy,  is  used  to  characterize  a  broad  spectrum  of 
steels  which  may  have  widely  varying  chemical  compositions  and  mechanical  properties. 
HSLA  steels  generally  have  higher  yield  strengths  than  plain  carbon  structural  steels  and 
may  also  show  varying  amounts  of  other  beneficial  properties  such  as:  greater  fracture 
toughness,  formability,  weldability,  and  good  atmospheric  corrosion  resistance  [Ref.  1]. 
Several  categories  of  HSLA  steels  include:  weathering  steels,  control-rolled  steels, 
pearlite-reduced  steels,  microalloyed  steels,  acicular  ferrite  steels,  and  dual  phase  steels 
[Refs.  1 ,  2:  p.  4*50  ]. 

Microalloying  in  the  HSLA  steels  is  capable  of  producing  many  characteristics  desired 
for  naval  shipbuilding  applications.  The  very  low  carbon  content  in  some  of  these  steels  in 
conjunction  with  niobium  reduces  the  formation  of  iron  carbides  and  instead  form  niobium 
carbonitrides  during  hot  rolling.  The  carbonitrides  retard  grain  growth,  providing  grain 
refinement,  increasing  strength  and  toughness  without  increasing  the  ductile-to-brittle 
transition  temperature  (DBTT)[Refs.  1 ,  3].  Molybdenum  increases  the  aged  strength  of  the 
steel,  improves  toughness  ,  and  hardenability  while  reducing  grain  boundary 
segregation[Ref.  4]. 

Copper  in  concentrations  of  greater  than  0.75  percent  increases  the  aged  strength  of 
HSLA  steels.  On  aging  at  temperatures  over  500°C  (930°F)  coherent  body  centered  cubic 
clusters  of  copper  form,  precipitation  hardening  the  steel.  On  continued  or  subsequent  aging 

these  coherent  particles  transform  into  non-coherent  e-phase  Cu  particles.  At  the  point  of 


maximum  strength  only  a  few  particles  are  visible  under  an  electron  microscope  [Refs.  3, 
5].  They  are  hard  to  "see"  because  they  are  coherent.  The  particle  size  at  maximum 
strength  has  been  shown  to  be  50  angstroms.  Copper  is  also  the  element  responsible  for 
increasing  the  steels  atmospheric  corrosion  resistance. 

Problems  associated  with  copper  as  an  alloying  element  are  related  to  its  relatively  low 
melting  temperature,  1083°C  (1981°F)  for  Cu  versus  1536°C  (2798°F)  for  Fe.  This 
becomes  a  problem  when  solution  heat  treating  a  copper  precipitated  steel  or  when 
hot-working  the  steel  at  a  temperatures  above  1050°C  (1922°F)  with  copper 
concentrations  above  0.5  percent.  Copper  builds  up  in  steel  scrap  and  cannot  be  removed  in 
subsequent  steel  production  heats  creating  further  difficulties  in  reprocessing.  Nickel  is 
added  to  the  steel  to  prevent  the  molten  copper  from  penetrating  the  austenite  grain 
boundaries  (hot-shortness).  [Refs.  3,  5] 

One  of  the  major  beneficial  qualities  of  HSLA  steels  is  its  weldability  without  any 
pre-heating,  due  primarily  to  the  low  carbon  content  of  the  steels.  Some  of  the  weld  metal  / 
heat  affected  zone  (HAZ)  characteristics  include  good  ductility  and  toughness  without 
significant  losses  in  strength  compared  to  the  base  metal.  [Refs.  1 ,  6] 

B.    DEVELOPMENT  OF  HSLA-1 00 

Weld  preparation  is  a  very  high  cost  and  time  consuming  process  in  naval  ship 
construction.  The  HY  series  of  steels  exhibit  both  the  required  strength  and  toughness  in 
base  metal  and  welds  but  entail  detailed  and  expensive  welding  procedures.  An  already 
qualified,  low  carbon,  copper  precipitation  strengthened,  HSLA-80  steel  showed  excellent 
weld  strength  and  toughness  without  any  preheat  in  plates  less  than  3/4  inch  thick.  To 
expand  on  the  capabilities  demonstrated  by  this  steel,  a  research  program  funded  under  the 
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SSN-21  project  was  initiated  by  the  Naval  Sea  Systems  Command  and  David  Taylor  Naval 
Ship  Research  and  Development  Center.  This  project  was  completed  under  contract  to  the 
AMAX  Materials  Research  Center  (AMRC),  subcontracted  to  Phoenix  Steel  Corporation, 
Claymont  Delaware.  [Ret.  7] 

The  purpose  of  the  development  project  was  to  modify  the  current  HSLA-80  steel  into  a 
highly  weldable,  copper  precipitation  strengthened  steel,  with  carbon  content  less  than  0.07 
percent ,  sulfur  less  than  0.008  percent,  phosphorous  less  than  0.010  percent,  containing 
manganese,  nickel,  chromium,  and  molybdenum.  In  addition  the  steel  was  required  to  be 
solution  treated,  quenched  and  aged  with  a  nominal  0.2  percent  offset  yield  strength  ranging 
between  100  ksi  (689  MPa)  and  120  ksi  (827  MPa),  minimum  average  transverse  Charpy 
V-notch  impact  properties  of  55  ft-lb  (75  J)  at  0°F  (-18°C)  and  30  ft-lb  (41  J)  at 
-120°F  (-84°C),  minimum  elongation  of  18  percent  in  2  inches  (51  mm),  and  minimum 
reduction  in  area  of  45  percent.  These  requirements  dictate  that  strengthening  from  solid 
solution,  grain  refinement,  and  precipitation  must  be  present.  [Ref.  7] 

Three  phases  of  development  for  the  new  steel  were  defined  by:  laboratory  development 
of  prospective  compositions  and  heat  treatments  (Phase  I),  production  and  testing  of  a  trial 
commercial  heat  (Phase  II),  and  procurement  of  a  larger  quantity  of  steel  (Phase  III)  based 
on  the  results  of  testing  of  the  Phase  II  steel.  The  resultant,  interim,  specifications  for  a 
HSLA-100  trial  heat  were  as  specified  in  Appendix  A.  [Ref.  7] 

The  microstructure  of  the  selected  steel  is  a  100  percent  low-carbon  banite  with  large 
amounts  of  copper  (1 .6  percent)  and  nickel  (3.5  percent)  for  strength  and  toughness. 
Minimal  quantities  of  niobium  are  for  austenite  grain  size  control.  [Ref.  7] 

Sensitivity  to  quenching  from  the  solution  treating  temperature  was  found  to  be  a 
concern  only  in  plates  greater  than  one  inch.  The  effect  of  a  slow  quench  possibly  caused  by 
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an  overheated  quenching  media  would  be  inadequate  retention  of  copper  in  solid  solution  and 
lower  yield  strengths.  The  100  percent  banite  microstructure  was  producible  with  even  the 
slowest  quenches  used  during  the  development.  [Ref.  7] 
C.   PURPOSE  AND  SCOPE  OF  PRESENT  WORK 

Brittle  fracture  is  type  of  failure  characterized  by  rapid  crack  propagation  and  little  or 
no  associated  plastic  deformation.  There  is  a  generally  high  tensile  stress  and  lower  energy 
associated  with  this  type  of  fracture  initiation  when  compared  to  ductile  fractures.  Ductile 
fracture  is  characterized  by  microvoid  initiation,  growth,  and  coalescence,  resulting 
typically  in  a  cup  and  cone  morphology.  Initiation  sites  have  been  associated  with  decohesion 
of  carbides  or  inclusions.  [Ref.  8] 

The  type  of  fracture,  brittle  or  ductile,  depends  on  many  variables  such  as  chemical 
composition  (alloying)  and  processing  of  the  steel.  Other  dependent  variables  include  the 
section  size,  strain  rate,  and  the  temperature  at  which  the  test  was  conducted.  Higher 
temperatures  are  generally  associated  with  ductile  fractures  and  lower  temperatures  with 
brittle  fractures.  There  is  a  temperature  range  associated  with  the  transition  from  ductile 
to  brittle  type  fractures,  denoted  by  a  temperature  called  the  ductile  to  brittle  transition 
temperature  (DBTT).  [Ref.  8] 

There  are  several  models  used  to  describe  crack  initiation  and  propagation  for  cleavage 
fractures  but  they  rely  heavily  upon  iron  carbides  as  crack  initiation  sites.  HSLA-100 
contains  very  low  levels  of  carbon  and  these  models  have  difficulty  describing  initiation  of 
cleavage  in  the  upper  transition  temperature  regions.  To  better  understand  and  model  the 
toughness  behavior  of  HSLA-100  and  future  steels  a  program  is  in  progress  to  investigate 
the  micro  mechanics  of  HSLA-100  fracture  behavior.  In  the  first  of  three  stages  of  this 
project  one  of  the  requirements  is  to  derive  true  stress  /  true  strain  tensile  curves  as  a 
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function  of  temperature  which  will  provide  the  necessary  Theological  input  for  a  finite 
element  analysis  of  HSLA-100.  The  purpose  of  this  research  to  derive  this  constitutive 
equation  of  the  stress  as  a  function  of  strain  and  temperature  under  the  conditions  of  a 
quasi-static  deformation.  [Ref.  8] 
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II.  BACKGROUND 


A.  TENSILE  TESTING 

The  data  obtained  in  a  tensile  test  is  generally  load,  measured  by  a  load  cell  or  some  other 
means,  and  displacement,  measured  by  an  extensometer  or  crosshead  movement.  When 
plotted  in  typical  fashion  with  load  plotted  on  the  ordinate  (Y-axis)  and  displacement  on  the 
abscissa  (X-axis)  the  general  form  of  the  curve  is  as  shown  in  Figure  1 . 


Max  Load 


Fracture 


Displacement 


> 


Figure  1. 
Typical  Load/Displacement  Curve 


In  this  Figure  there  are  three  important  points.  First  is  the  yield  point  where  the  curve 
deviates  from  the  initial  approximately  linear  portion  of  the  curve.  The  point  of  maximum 
load  is  the  second,  where  necking  begins  and  corrections  for  the  triaxial  stress  state  at  the 
minimum  diameter  begins.  The  curve  from  this  point  shows  a  decreasing  load  carrying 
ability  as  displacement  continues,  until  fracture  when  the  specimen  incapable  of  carrying 
any  load.  [Ref.  9:  pp.  13-15] 
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Conversion  of  the  load  into  true  stress  and  displacement  into  true  strain  results  in  a 
more  usable  form  of  the  information  for  engineering  purposes.  True  stress  is  calculated  by- 
the  equation: 

aT=  F/A; 

where  Fj  is  the  instantaneous  load  carried  by  the  specimen,  and  Aj  is  the  instantaneous 

cross-sectional  area  of  the  specimen.  The  true  stress,  Oj,  describes  the  state  of  stress,  in 

units  of  force/unit  area,  in  the  plane  where  that  cross  sectional  area  was  measured. 

A  fundamental  assumption  associated  with  calculating  the  true  strain  is  that  the  volume 
of  the  sample  remains  constant  throughout  the  deformation.  With  this,  volume,  V= 

constants  A-,  lj=  AQ  lQ,  where  lj  is  the  instantaneous  length  and  the  subscript  'o'  denotes  the 

initial  conditions.  True  strain  is  defined  as: 

€T=  Ln(  VAj)-  Ln(lj/I0)=  2  Ln(d0/dj) 

A  typical  true  stress-true  strain  curve  is  diagrammed  in  Figure  2. 

In  this  curve  the  section  from  the  origin  to  point  A  is  the  elastic  portion  of  the  curve. 
The  elastic  portion  of  the  stress-strain  curve  is  nearly  linear  and  thus  the  stress  strain 
ratio  should  be  a  constant.  This  ratio  was  initially  described  by  Robert  Hooke  in  1678  and 
bears  his  name,  Hooke's  law  [Ref.  10:  p.  6]. 

E=  a/e . 
In  this  equation  the  constant  E,  called  Young's  modulus,  is  one  of  the  fundamental  constants 
describing  the  stiffness  of  the  material.  To  define  the  yield  point,  frequently  a  difficult 
point  to  specify  exactly,  another  commonly  used  property  is  defined  as  the  0.2  percent  offset 
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yield  strength.  To  arrive  at  the  0.2  percent  offset  yield  strength  a  line  is  drawn  parallel  to 
the  elastic  portion  of  the  stress-strain  curve,  slope  equal  to  E,  intersecting  the  X-axis  at 
0.2  percent  strain.  Where  the  line  crosses  the  stress  /  strain  curve  defines  the  0.2  percent 
offset  yield  strength. 


•K 


True 
Stress 


Bridgman  Corrected  True  Stress 


» 


True  Strain     € 
Figure  2. 
Typical  True  Stress-True  Strain  Curve 

The  second  portion  of  the  curve,  from  point  A  to  B,  is  the  region  of  uniform  plastic 

deformation  and  generally  follows  an  equation  of  the  form: 

€n=  a/K. 
The  final  portion  of  the  stress-strain  curve,  from  point  B  to  failure,  is  called  the  region  of 
non-uniform  plastic  deformation  and  once  again  strain  is  approximately  proportional  to 
stress.  The  lower  curve  shown  in  Figure  2  describes  the  stress  corrected  for  the  triaxial 
stress  state  which  exists  in  the  neck.  This  will  be  discussed  in  more  detail  later.  [Ref.  9:  pp 
13-17] 
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In  the  latter  two  portions  of  the  curve  there  are  two  strain  hardening  regimes.  Strain 
hardening  is  the  phenomenon  whereby  a  material  strained  past  the  yield  point,  when  relaxed 
will  follow  a  line  approximately  parallel  to  the  elastic  portion  of  the  curve  and  when 
retensioned  will  begin  to  yield  at  a  load  higher  than  the  initial  yield  point.  The  difference  in 
the  yield  points  is  the  strain  hardening  and  the  amount  of  strain  not  recovered  during 
relaxation  is  the  plastic  deformation.  Expressing  mathematically  the  plastic  deformation  as 


a  strain,  the  plastic  strain,  e  ,  is 


€p=et-(a/E) 

where  the  quantity  a/E  is  the  elastic  strain  of  the  material.  [Ref.  9:  pp.  15-17] 

When  these  stress-strain  points  are  plotted  on  logarithmic  coordinates  the  shape  of  the 
curve  will  often  be  a  shallow  S-shape  as  shown  in  Figure  3. 
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Log 

True 
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Log  True  Strain 

Figure  3. 
Typical  Logarithmic  Stress-Strain  Curve 

This  form  of  the  stress-strain  curve  is  the  basis  of  two  stress-strain  relationships  which 

will  be  discussed  in  more  detail. 
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B.   STRESS  STRAIN  RELATIONSHIPS 

Over  the  years  many  forms  of  stress-strain  relations  have  been  formulated.  Two  which 
have  gained  some  recognition  are  the  Power  Equation  and  the  Voce  equation. 

1.    The  Power  Equation 

The  power  equation,  attributed  to  Ludwik  (1909)  (a=  a    +  Ken)  and  Hollomon 

(1945)  (g=  Ken),  were  widely  used  for  many  years  to  describe  stress-strain  behavior. 
When  the  logarithmic  plot  of  the  stress-strain  curve  is  viewed  on  some  scales,  a  straight 
line  appears  to  fit  the  curve  very  well  through  all  portions  of  the  curve.  The  Hollomon 
version  of  the  equation  is  more  widely  used  and  may  take  either  of  the  forms: 

o=  Kem     or: 

log  o=  log  K  +  m  log  €. 
K  is  called  the  strength  coefficient  and  represents  the  true  stress  at  a  unit  true  strain.  The 
quantity  log  K  is  the  least  squares  Y-intercept  of  the  straight  line  representation  of  the 
curve.  The  slope  of  the  line,  m,  is  the  slope  of  the  line  on  logarithmic  coordinates  and 
represents  the  strain  hardening  rate.  It  is  mathematically  equal  to  the  true  strain  at  the 
necking  point  (maximum  load).  The  coefficients  for  this  equation  are  affected  by  variables 
of  testing  including  temperature  and  strain  rate.  [Ref.  1 1 :  pp  155-156] 

This  equation  has  been  shown  to  have  limited  success  in  describing  stress-strain  data 
beyond  the  point  of  necking  (maximum  load).  This  equation  is  more  successful  when  used 
with  true  plastic  strain  instead  of  the  true  strain.  Because  of  the  simplicity  of  the  power 
equation  it  is  good  for  fast  comparisons  of  materials.    The  solution  of  the  coefficients  can  be 
easily  determined  by  a  least  squares  linear  curve  fit  of  logarithmic  stress-strain  data 
points.  The  slope  of  the  line  is  the  strain  hardening  coefficient  and  the  intercept  is  equal  to 
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the  log  of  the  strength  coefficient.  [Refs.  1 0:  pp  373-375, 1 1 :  pp1 55-1 56] 

2.    The  Voce  equation 

Voce  disagreed  with  the  form  of  the  power  equation  because  of  its  empirical  nature 
and  its  failure  to  accurately  describe  the  S-shaped  form  of  logarithmic  stress-strain  curve. 
In  addition,  the  power  equation  implies  that  load  carrying  ability  or  strain  hardening  would 
infinitely  increase  as  the  material  is  deformed.  What  is  seen  experimentally  is  a  leveling 
off  at  some  value  for  incremental  increases  in  strain.  [Ref.  1 1 :  pp.1 60-1 69, 1 2] 

To  arrive  at  a  different  form  of  an  equation  describing  the  stress-strain  curve  Voce 

iooked  at  the  logarithmic  curve  for  copper  samples  and  picked  the  approximately  asymptotic 

values  at  small  and  large  strains.  Several  terms  Voce  used  in  his  equation  and  some 

alternate  forms  of  them  follow: 

TABLE  1      SYMBOLOGY  USED  WITH  THE  VOCE  EQUATION 
Primary  Alternate  Definition 

Sb  o~o  The  initial  or  threshold  stress  where  homogeneous  plastic 

deformation  begins. 
S°o  aoo  The  final  constant  stress  at  which  deformation  proceeds  without 

an  increase  in  stress. 
Cb  Equal  to  (S°o  -  So),  the  plastic  stress  capacity  of  the  material. 

C  Equal  to  (S<»  -S),  the  plastic  stress  remaining  in  the  material. 

R  The  strain  ratio,  equal  to  l/lo  =  Ao/A  =  (do/d)*  in  tension 

testing. 

Ln  R=  €. 
A  €_,  k  A  constant  called  the  characteristic  strain. 

In  his  analysis  of  copper  samples  in  compression  he  noted  that  a  plot  of  the  strain  versus 
Ln(C/Co)  resulted  in  a  linear  plot.  By  rearranging  the  terms  and  solving  for  stress  he 
arrived  at  the  final  form  of  the  Voce  equation: 

S=  Soo  -(Soo  -  So)  e'€/A 
This  form  of  the  equation  has  equal  or  better  results  in  fitting  stress-strain  data  than  the 
power  equation  and  in  most  cases  is  applicable  for  tensile  tests  as  well  as  compressive  tests. 
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It  can  be  used  for  any  material  showing  the  characteristic  curve  shown  in  Figure  3.  A 

problem  associated  with  the  use  of  this  equation  is  that  it  is  most  accurate  at  high  strains 

thus  does  not  accurately  describe  the  elastic  portion  of  the  stress-strain  curve  and  may  fit 

poorly  in  the  region  of  uniform  plastic  deformation.  [Refs.  1 1 :  160-162, 12] 

The  Voce  model  does  have  some  physical  reality  (thus  not  purely  empirical).  For 

example,  if  there  is  some  finite  number  of  strain  hardening  sites  and  a  small  increment  of 

strain  is  applied  to  the  material,  a  small  portion  of  the  strain  hardening  sites  will  be 

blocked.  This  proportion  of  the  blocked  sites  to  the  total  number  of  hardening  regions  is 

equal  to  the  proportion  of  the  plastic  stress  capacity  used  in  the  material  as  compared  to  the 

total  plastic  stress  capacity  of  the  material: 

n   _  (S  -  So^ 
N       (Soo-S) 

Viewed  from  a  dislocation  movement  point  of  view,  the  initial  asymptotic  portion  of  the 

curve  exists  because  dislocations  are  free  to  move  with  none,  or  very  few  being  blocked.  As 

deformation  proceeds  an  increasing  amount  of  the  dislocations  become  blocked  until  finally 

the  stress  asymptotically  approaches  S°°  because  the  blockage  rate  and  annihilation  rate  of 

dislocations  balance.  [Refs.  11:  pp.  164-165, 12] 

In  the  treatment  of  specimens  which  have  some  unknown  degree  of  plastic  deformation  in 

them,  the  power  equation  is  incapable  of  detecting  the  initial  strain  hardening  associated 

with  it.  If  the  Voce  equation  coefficients  are  known  for  the  material,  the  amount  of  prestrain 

on  the  material  can  be  predicted  accurately  from  the  point  where  the  prestrained  specimen 

meets  the  original  Voce  curve,  The  value  of  A,  and  S°o  will  remain  unchanged.  [Ref  1 1 :  p. 

171] 
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C.   CORRECTION  FOR  THE  TRIAXIAL  STRESS  STATE  IN  THE  NECK 

The  maximum  load  point  of  the  load-displacement  curve  is  the  point  where  necking 
begins.  The  axial  stress  state  in  the  neck  and  the  load  measured  by  the  load  cell  are  not  equal 
because  of  a  radial  component  of  flow  for  material  in  the  vicinity  of  the  neck  of  a  cylindrical 
sample.  The  physicist,  Percy  Williams  Bridgman,  studied  on  this  problem.  In  his  analysis 
of  the  true  axial  stress  in  the  neck  of  a  tensile  specimen  he  used  tensors  in  cylindrical 
coordinates,  equilibrium  and  boundary  conditions  (compatibility)  at  the  surface  as  a  start 
point  in  the  analysis.  The  round  tensile  samples  are  assumed  to  have  rotational  symmetry 
and  picking  the  origin  at  the  minimum  diameter  or  the  sample,  there  is  assumed  to  be 
symmetry  on  the  ±  z-axis.  Applying  the  von  Mises  criterion  of  maximum  distortion  energy 
to  the  stresses  he  was  able  to  arrive  at  an  equation  which  approximates  the  actual  uniaxial 
stress  state  in  the  neck. 

<?=  aa„  /  (1  +  4R/d)  Ln(1  +  d/4R) 

dV 

In  this  equation,  a  is  the  stress  corrected  for  the  triaxial  stresses  in  the  neck  and  is  always 
less  than  the  stress  measured  by  the  load  cell  (divided  by  the  area).  The  measured  stress  is 

g_v,  R  is  the  radius  of  curvature  in  the  vicinity  of  the  minimum  diameter  of  the  neck,  and  d 
is  the  minimum  diameter  of  the  neck.  The  quantity  ,  1  /  (1  +  4R/d)  Ln(1  +  d/4R),  is 

equal  to  a  /  a  av  and  is  the  Bridgman  correction  [Refs.  10:  pp.  566-569, 13]. 

In  an  attempt  to  generate  an  empirical  relation  of  the  Bridgman  correction  to  more 
readily  obtained  stress-strain  information  Bridgman  and  others  had  some  success  in 
correlating  the  Bridgman  correction  to  the  true  strain  minus  the  true  strain  at  the  necking 

point,  (e  -  €  n).  The  curves  they  generated  were  satisfactory  for  approximate  corrections 
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but  in  order  to  determine  the  actual  correction  the  actual  ratio  of  d/R  must  be  measured. 
[Refs.  9:  pp.  21-22, 13] 

D.  TEMPERATURE  EFFECTS  ON  MATERIAL  BEHAVIOR 

Strength  and  plastic  deformation  are  temperature  dependent  phenomenon  with  the 
controlling  mechanisms  existing  at  the  crystallographic  scale.  Resistance  to  plastic 
deformation  is  determined  by  the  shear  stress  required  to  make  dislocations  glide  or  move  in 
their  slip  planes.  There  are  temperature  and  strain  rate  dependent  components  of  this 
stress.  The  stress  applied  is  required  to  enable  the  dislocation  to  overcome  barriers  to 
movement  such  as  precipitates,  grain  boundaries,  and  other  dislocations.  These  examples 
are  considered  long  range  fields.  Another  basic  type  of  barrier  to  movement  in  any  crystal 
lattice  is  a  short  range  barrier  called  a  Peierls-Nabarro  barrier.  This  barrier  results 
from  the  inherent  resistance  of  the  crystal  lattice  to  movement  of  dislocations.  [Ref.  10: 
320-326] 

A  dislocation  is  a  defect  in  a  crystal  lattice  which  arises  from  the  addition  or  omission  of 
a  planar  series  of  atoms.  Movement  of  a  dislocation  comes  about  from  a  movement  of  atoms 
in  the  lattice.  The  Peierls-Nabarro  barrier  can  be  thought  of  in  a  simplistic  manner  as  the 
resistance  offered  by  neighboring  atoms  to  the  movement  of  the  atom  in  an  ordered 
structure.  The  energy  required  to  overcome  this  barrier  is  called  the  Peierls-Nabarro 
energy.  At  temperatures  near  room  temperature  (20°C)  this  energy  is  predicted  to  be 
relatively  low  but  at  low  temperatures,  in  materials  with  strong  directional  bonding 
including  iron,  these  stresses  increase  exponentially  as  absolute  zero  (0°K)  is  approached 
and  are  responsible  for  a  significant  portion  of  the  yield  strength  of  the  material.  [Ref.  10: 
pp.  324] 
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III.  EXPERIMENTAL  PROCEDURE 

A.   MATERIAL  AND  SAMPLE  PREPARATION 

A  section  of  1-1/4  inch  (32  mm)  HSLA-100  plate  (Plate  #  5644-1 6B)  from  the  trial 

commercial  heat  of  HSLA-100  was  provided  by  the  David  Taylor  Naval  Ship  Research  and 

Development  Center  for  use  in  this  research.  The  plate  was  heat  treated  by  the  supplier  by 

solution  treating  at  949°C  (1650°F)  for  70  minutes  and  water  quenching  and  then  aging  at 

566°C  (1050°F)  for  70  minutes  and  water  quenching.  The  strength  properties  listed  by 

the  supplier  are  listed  in  Table  2. 

TABLE  2 

STRENGTH  PROPERTIES  OF  HSLA-100  PLATE  #5644-1 6B 

(AS  REPORTED  BY  THE  SUPPLIER) 


Yield 
Strength 
(ksi)  " 

Ultimate 
Tensile 
Strength 
(ksi) 

Percent 
Elongation 

Percent 
Reduction 
in  Area 

Top 
Transverse 

101 

147 

22 

65 

Bottom 
Transverse 

106 

139 

23 

65 

Reference  14  noted  the  material  properties  observed  in  his  research  deviated  significantly 
from  those  listed  in  Table  2  with  much  higher  yield  and  tensile  strengths.  This  was 
confirmed  separately  by  David  Taylor  NSRDC  and  through  subsequent  testing  of  various 
additional  heat  treatments  a  determination  was  made  to  additionally  heat  treat  this  steel  at 
621  °C  (1150°F)  to  bring  the  yield  strength  within  the  desired  range  of  100  to  120  ksi 
(689  to  758  MPa). 
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Tensile  blanks,  measuring  approximately  5/8  inch  square  by  4-1/2  inch  long,  were  cut 
from  the  plate  transverse  to  the  rolling  direction.  Three  tensile  blanks  were  wired  together 
with  spacers  and  an  additional  blank  for  internal  temperature  monitoring.   The  set  of 
tensile  blanks  was  then  placed  in  a  furnace  stabilized  at  621  °C  (1150°F)  and  aged  at  621  °C 
±  1  percent  for  60  minutes  and  water  quenched.  The  oxides  that  formed  on  the  surface  of  the 
bar  during  the  aging  process  were  removed  from  the  face  of  the  bar  perpendicular  to  the 
rolling  direction,  and  the  hardness  was  determined  on  the  Rockwell  C  scale  by  the  average  of 
5  measurements.  The  average  of  all  52  blanks  aged  was  RC  24.8  with  a  standard  deviation 
of  0.98,  corresponding  to  the  accuracy  of  the  testing  apparatus  used.  This  average  hardness 
converts  to  an  approximate  tensile  strength  of  118  to  122  ksi  (814  to  841  MPa).  The 
resultant  microstructures  are  shown  in  Figures  4  and  5. 

In  the  light  micrograph  of  Figure  4,  the  microstructure  is  seen  as  a  fine  banite  with 
banding  lightly  evident  from  this  aspect,  looking  at  the  side  of  the  plate,  an  end  view  of  the 
samples.  Figure  5  is  a  TEM  micrograph  illustrating  high  dislocation  density  and  copper 
precipitates  evident  within  grains  and  at  the  grain  boundaries,  visible  as  oblong  spots. 

On  completion  of  the  aging  two  blanks  were  machined  into  uniform  gage  length  samples  in 
accordance  with  Figure  6  (a)  and  the  remainder  of  the  tensile  blanks  were  machined  into 
hourglass  shaped  specimens  in  accordance  with  Figure  6  (b).  The  diameters  of  the  finished 
tensile  specimens  were  then  measured  and  recorded. 
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Figure  4 
Light  Micrograph  of  Aged  HSLA-100  •  2  percent  Nital  Etch  for  25  seconds 

500X 


Figure  5 
TEM  Micrograph  of  Aged  HSLA-100 
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B.   TEST  APPARATUS 

Tensile  tests  were  conducted  with  a  MTS  810  Material  Test  System  including  as  the 
major  functional  parts  (all  manufactured  by  MTS):  409  temperature  controller,  464  data 
display,  410  digital  function  generator,  445  controller,  413  master  control  panel, 
Thermotron  environmental  control  chamber,  10  ton  load  cell,  632.1 9B  diametral 
extensometer,  and  load  frame  with  hydraulic  actuators.  The  diametral  extensometer  was 
required  to  follow  the  deformation  in  the  neck  to  failure  using  the  modified,  smaller  radius, 
tips  described  in  Reference  14. 

Several  modifications  to  the  basic  to  the  basic  setup  were  made  to  improve  temperature 
control.  Phenolic  washers  were  used  in  all  unions  in  the  actuator  train  to  restrict  heat  flow 
into  the  system.  For  cooling  to  the  sample  grips,  the  liquid  nitrogen  tubing  was  split  so  that 
both  grips  would  receive  equal  flow  of  nitrogen  and  the  tubing  size  was  increased  to  3/8  inch 
from  1/4  inch  to  allow  the  nitrogen  gas  which  boiled  off  to  escape  easier.  Without  the 
larger  tubing  the  backpressure  created  by  the  nitrogen  boiling  off  could  shut  off  the  flow  of 
nitrogen  to  the  grips.  These  modifications  resulted  in  slower  temperature  changes  and 
allowed  faster  response  to  the  changes,  thus  better  temperature  control.  With  this  setup 
alone  temperatures  to  -180°C  (-308°F)  were  easily  attained  and  controlled.  Photographs 
of  the  MTS  control  system,  grips  and  chamber  are  shown  in  Figure  7  (a)  and  (b). 
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NOTES:  1.  All  dimensions  in  inches. 

2.    Tolerances  as  per  ASTM  tensile  specimen  standards. 

3.  Specimen  gage  length  to  be  perpendicular  to  plate  as  rolled  direction. 
4.  Gage  length  shall  be  32  rms. 
5.  Mark  with  applicable  specimen  number  on  both  ends. 


5  in  |.5  in    .5  in 
* — ► 


4.00  in 
* — 1.0  ■ 


I — 1/2-20    NF 


.25 
2.75  R 


(a.) 


NOTES:  1.  All  dimensions  in  inches. 

2.  Tolerances:  As  per  ASTM  tensile  specimen  standards. 
3.  Specimen  gage  length  to  be  perpendicular  to  plate  as  rolled  direction. 
4.  Reduced  section  area  of  specimen  shall  be  polished  in  a  manner 
parallel  to  specimen  longitudinal  axis  to  32  rms. 
5.   Mark  with  applicable  specimen  number  on  both  ends,  vibrating  type 
engraving  tool  is  permissible. 


4.00 


1.02  — ►[ 


1.1  R 


1/2-20    NF 


(b.) 

Figure  6 
Tensile  Specimen  Drawings 
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Figure  7 
MTS  Function  Generator,  Controls  and  Test  Assembly 

Temperature  control  was  achieved  by  two  different  methods  dependent  on  the 
temperature  desired.  At  temperatures  below  -150°C  (-238°F)  the  cooling  circuit  switch 
had  to  be  manually  cycled  and  at  temperatures  above  that  the  installed  temperature 
controller  was  capable  of  maintaining  the  temperature  within  ±2°C.  The  top  grip 
temperature  was  monitored  internally  and  was  the  sensing  point  for  the  controller. 
Temperature  was  sensed  by  five  chromel-alumel  thermocouples  sensing  chamber  air 
temperature,  top  and  bottom  grip  temperature,  and  sample  temperature  above  and  below  the 
hourglass  gage  length.  The  thermocouples  were  calibrated  in  liquid  nitrogen  indicating 
-189°C  for  an  actual  temperature  of  -1 95.8°C. 

Procedures  for  use  of  the  tensile  test  system  and  general  conduct  of  testing  are  those 
provided  by  MTS  for  this  system  modified  for  the  specific  requirements  of  these  tests  by 
Hamilton  in  Reference  14.  An  additional  check  sheet  was  used  to  emphasize  particular  parts 
of  the  testing  process. 
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Figure  8 
Environmental  Test  Chamber,  Actuator,  and  Grips 

C.   DATA  COLLECTION  AND  REDUCTION 

During  testing,  progress  of  the  test  was  monitored  and  recorded  using  an  X-Y  plotter 

sensing  the  displacement  and  load  voltages.  This  plot  was  an  important  indication  of 

experimental  irregularities  and  was  used  to  verify  information  gathered  by  other  means. 

The  load  and  displacement  voltages  were  sensed,  converted,  and  recorded  using  a 

Hewlett-Packard  3437A  system  voltmeter,  and  3497  data  acquisition/control  unit,  both 

controlled  by  a  Hewlett-Packard  9826  computer  all  shown  in  Figure  9.  Software  for  the 

data  acquisition  process  was  a  truncated  form  of  'JHCOLLECT"  written  by  Hamilton  and 

provided  in  Reference  14.  Modifications  to  the  program  included  sampling  only  load  and 
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displacement  voltages,  increasing  the  size  of  the  arrays  for  data  storage,  and  sampling  at  a 
constant  frequency  throughout  the  test.  The  increased  size  of  arrays  was  required  to  store 
large  amounts  of  data  collected  at  intervals  of  less  than  a  second.  The  uniform  sampling 
provided  that  no  one  portion  of  the  curve  would  be  weighted  differently  in  a  least  squares 
reduction  routine. 

The  large  data  files,  containing  approximately  1 800  to  2200  data  sets,  were  stored  on 
disc  in  the  Hewlett-Packard  computer  then  retrieved  at  some  linear  interval  to  select 
approximately  1 1 0  load/displacement  data  points  from  the  start  of  the  test  until  failure. 
These  data  points  were  then  manually  transferred  into  a  Macintosh  personal  computer  for 
analysis  and  plotting  of  the  resultant  curves. 

Fracture  surfaces  were  examined  by  scanning  electron  and  light  microscopy  for 
examination  of  failure  types.  Metallographic  examination  was  also  conducted  by  light 
microscopy. 


Figure  9 
Data  Acquisition  System 
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IV.  ANALYSIS  AND  RESULTS 

A.  DETERMINATION  OF  THE  ELASTIC  MODULUS 

In  order  to  convert  true  strain  to  true  plastic  strain,  the  elastic  modulus  for  the 
material  was  first  determined.  The  sole  purpose  for  the  uniform  gage  length  samples  was  to 
achieve  this  end.  The  Young's  modulus  for  the  re-aged  HSLA-100  in  the  transverse 
direction  was  determined  using  the  two  uniform  gage  length  samples  machined  from  the  plate 
in  the  transverse  orientation,  longitudinal  extensometer  and  the  experimental  setup  as 
described  previously.  Each  of  the  two  samples  were  loaded  slowly  to  what  was  expected  to  be 
about  60  percent  of  the  yield  strength,  3.2  kips.  The  voltages  were  recorded  on  the  X-Y 

plotter,  data  converted  to  stress/strain  and  the  slope  for  each  of  the  two  samples,  da/de, 
determined.  For  the  two  samples,  specimen  numbers  16  and  17,  the  elastic  modulus  were 
found  to  be  24,  053  ksi  and  23,  953  ksi  respectively  for  an  average  Young's  modulus  of 

E=  24,000  ksi  (1654.7  MPa). 
This  value  is  close  to  that  reported  by  Hamilton  for  the  as  received  HSLA-100  (  24,140 
ksi)[Ref.  14]. 

Because  of  the  expanded  scale  of  the  plot  required  for  the  determination  of  the  the  elastic 
modulus  improved  chamber  supports  were  manufactured  to  remove  a  small  discontinuity  in 
the  plot  at  very  low  loads  and  strains. 

B.  INITIAL  TESTING  AND  ANALYSIS  OF  HOURGLASS  SPECIMENS 

Four  samples,  numbers  39  to  42,  were  initially  tested  to  the  practical  full  travel  of  the 
diametral  extensometer,  about  .077  inches  (  0.736  in/in  true  strain),  in  strain  control 
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then  manually  broken  in  stroke  control.  The  computer  data  acquisition  only  collected  the 
initial  full  range  travel  of  the  extensometer.  The  final  radius  of  curvature  of  the  neck  and 
final  diameter  were  measured  and  recorded.  The  data,  converted  to  the  form  of  load  and 
displacement,  was  stored  and  a  reduced  number  of  data  points  transferred  for  analysis. 
1.    Computer  Aided  Determination  of  the  Voce  Coefficients 

To  determine  the  coefficients  for  this  equation  is  difficult  since  three  unknowns  are 
involved.  However  using  one  of  the  derivative  forms  of  the  Voce  equation  results  in  an  easier 
solution.  Voce  noted  a  linear  relationship  when  he  plotted  the  true  strain  versus  the 
quantity  Ln[(S°o-S)/(S°°-So)].  The  characteristic  strain,  A,  would  be  the  slope  of  this 

line.  Since  the  intercept,  Ln(Soo-So),  is  a  constant,  a  plot  of  the  equation:  e=-A  Ln[S°°- 
S]  should  yield  a  straight  line,  but  this  equation  still  has  two  unknowns.  [Ref.  1 1 :  p.  161] 

The  recommended  solution  of  this  problem  in  is  to  select  a  true  stress  higher  than 
any  observed,  set  it  equal  to  Soo,  and  iterate  this  value  of  S°°  until  the  straightest  line 
value  for  Soo  is  found  using  a  comparative  method  such  as  the  correlation  coefficient  [Ref. 
11 :  pp.  1 66-1 69].  A  computer  program  in  Microsoft  BASIC  was  written  to  do  this. 
Several  versions  for  the  program  were  written  for  calculation  of  the  Voce  coefficients  for 
true  stress/true  strain,  Bridgman  corrected  true  stress/true  plastic  strain  and  so  on  but 
only  the  true  stress/true  strain  version  is  included  as  Appendix  B. 

Figure  10  shows  the  effect  of  changing  S°o  on  the  slope  (-A)  of  this  relationship.  As 
the  true  stress  approaches  Soo,  the  difference  between  Soo  and  the  true  stress,  S,  becomes 
smaller.  The  Ln  function  becomes  more  sensitive  to  fluctuations  in  the  load  resulting  in  the 
spikes  in  this  plot  which  are  amplified  at  higher  stresses.  The  gap  in  the  data  points  results 
from  stopping  and  restarting  the  test.  The  values  for  the  constants  derived  by  this  method 
are  very  sensitive  to  the  data  from  each  extreme  (high  and  low  stress)  of  the  plot.  The  data 
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below  the  yield  point  approaches  the  horizontal  so  inclusion  of  these  points  in  a  least  squares 
fit  will  pull  the  straight  line  fit  to  the  right  at  the  bottom,  lowering  the  value  of  the 
characteristic  strain  and  also  produce  a  lower  value  of  S<»  and  So  for  the  best  least  squares 
fit.  From  the  actual  data  set  of  Figure  1 0,  (sample  50)  the  final  value  of  S°o  was 
determined  to  be  265.0807  ksi,  a  characteristic  strain  of  0.7891093;  the  correlation 
coefficient  for  the  equation  using  these  values  for  the  contants  is  0.9968  (a  correlation 
coefficient  of  1 .0  indicates  a  perfect  linear  fit). 
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Figure  10 
Determination  of  S°°  and  the  Characteristic  Strain 


Analyzing  the  data  of  these  first  four  tests  (samples  39  to  42)  brought  about  several 
important  results.  The  first  problem  was  that  of  determining  the  actual  load/displacement 
at  failure.  To  determine  this  point  the  final  diameter  was  measured  from  the  failed  sample 
and  the  load  was  estimated  from  the  X-Y  plot.  The  results  of  sample  number  39  are  typical 
of  the  results  and  the  plot  of  true  stress/true  strain  is  shown  in  Figure  1 1 .  The  larger 
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points  correspond  to  the  first  point  considered  in  the  Voce  analysis  (approximately  the  yield 
point),  the  point  of  maximum  load  and  the  last  data  point  acquired.  The  point  on  the  top  right 
is  the  estimated  failure  point;  clearly  this  point  is  not  correct  because  an  extrapolation  of 
the  stress-strain  data  from  the  last  data  collected  would  not  even  come  close  to  this  point.  In 
addition  the  value  for  S°°,  where  the  true  stress  should  reach  a  final  asymptotic  value,  was 
calculated  as  226.4978  ksi,  far  below  the  maximum  true  stress  estimated  from  the 
fracture  point.  This  was  also  determined  to  be  low  because  the  lack  of  data  in  the  high  stress 
region  biases  the  linear  regression  to  the  low  stress  region  resulting  in  low  values  of  A  and 


HSLA-100  sample  39«  20°C 


300 


I  '  ■  ■  ■  I  ■  ■  ■  '  I 

0.0      0.1      0.2      0.3      0.4     0.5      0.6      0.7      0.8     0.9      1.0      1.1      1.2      1.3 

True  Strain  (in/in) 


Figure  1 1 
Estimate  of  the  Voce  Equation  and  Failure  point  -  Sample  Number  36 
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2.    Modifications  to  the  Testing  Procedure 

Based  on  these  results  it  was  determined  that  it  was  rrecessary  to  collect  strain  data 
right  up  to  fracture.  This  would  extend  the  curve,  giving  a  better  indication  of  the  actual 
stress  at  failure  and  properly  weight  the  least  squares  curve  fit  to  find  S°o. 

Another  series  of  measurements  and  analyses  being  conducted  at  the  same  time  dealt 
with  the  Bridgman  correction,  the  measurement  of  the  radius  of  curvature  and  how  it 
changes  with  the  strain.  Reference  9  [p.  21]  indicated  the  change  in  the  Bridgman 
correction  could  be  roughly  correlated  with  the  difference  between  the  instantaneous  true 
strain  and  the  strain  at  the  necking  point  (maximum  load)  and  that  this  could  be  used  for 
determining  approximate  Bridgman  corrections  as  a  function  of  strain.  To  confirm  this,  one 
test  was  conducted  where  the  loading  was  interrupted  at  several  points  between  the 
maximum  load  and  failure.  The  radius  of  curvature  and  diameter  would  be  measured  at  each 
interruption.  For  future  tests  it  was  also  decided  to  use  the  full  range  travel  of  the 
extensometer  (0.077  inches)  as  a  stop  point  for  measurement  of  the  radius  of  curvature 
giving  another  indication  of  the  trend  of  the  Bridgman  correction. 

To  accomplish  testing  in  this  manner  a  thin  (about  0.040  inch)  piece  of  copper  was 
shaped  to  fit  snugly  on  one  of  the  extensometer  tips.  This  additional  0.040  inch  would  be 
sufficient  to  extend  the  range  of  displacement  of  the  extensometer  to  accommodate  the  largest 
deformations  experienced  to  that  time.  To  correct  for  the  measured  displacement  change  that 
occurred  when  the  copper  foil  was  inserted,  this  displacement  (about  .040  inches)  was 
added  to  all  of  the  data  collected  during  the  retrieval  and  transfer  of  the  data  files. 

C.  COMPLETION  OF  TESTING  HOURGLASS  SPECIMENS 

A  total  of  22  additional  hourglass  shaped  specimens  were  tested  at  temperatures  varying 
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from  -1 89°C  to  37°C  (-308.2°F  to  98.6°F)  with  no  appreciable  difficulties  encountered. 
Temperature  was  maintained  at  ±5°C  at  temperatures  from  -150°C  to  -1 80°C  and  at  than 
±2°C  otherwise  and  the  strain  rate  was  constant  at  about  1x10"3  in/in  sec. 
1.    Voce  Coefficients  and  Comparison  to  the  Hollomon  Power  Equation 

The  Voce  equation  coefficients  were  determined  for  21  tests  using  the  computer 
program  described  in  Appendix  B.  A  typical  plot  of  the  load/displacement  is  shown  in  Figure 
12  (a.).  The  most  notable  feature  of  this  plot  is  the  manner  in  which  the  material  fails.  A 
sudden  drop  in  load  (at  a  displacement  of  0.1  inches)  occurred  just  before  the  final  fracture. 
This  phenomenon  was  noted  in  some  form  at  all  temperatures  greater  than  -100°C.  The 
true  stress/true  strain  plots  of  these  data  and  logarithmic  plots  are  shown  as  Figure  12  (b.) 
and  (a).  The  rapid  dropoff  of  the  load  is  also  seen  in  these  plots.  Another  feature  of  the 
logarithmic  plot,  which  is  slightly  visible,  is  a  small  discontinuity,  or  small  change  in  slope 
at  about  the  maximum  load  point.  Since  data  prior  to  the  yield  point  and  in  the  last  region  of 
the  rapid  drop  in  load  cannot  be  described  by  the  Voce  equation  they  were  not  considered  in 
the  Voce  analysis.  The  load/displacement  or  stress/strain  plot  had  to  be  viewed  first  to 
determine  the  region  which  could  be  described  by  the  Voce  Equation.  The  regions  not 
included  in  the  analysis  are  relatively  small  at  either  end  of  the  curves  and  more 
information  can  be  gained  by  their  exclusion. 
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Figure  12 
Typical  Experimental  Load/Displacement  and  True  Stress/True  Strain  Plots 
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Figure  12 
Typical  Experimental  Load/Displacement  and  True  Stress/True  Strain  Plots 


If  the  Hollomon  power  equation  is  appropriate,  the  log  true  stress  versus  log  true 
strain  should  be  a  straight  line.  When  a  plot  similar  to  Figure  12  (c.)  is  viewed  on  an 
expanded  scale  of  log  true  stress  ranging  from  zero  to  2.3  on  the  Y-axis  the  resultant  plot 
does  appear  to  be  quite  linear,  however  by  examining  the  plot  on  the  magnified  scale  as 
shown  in  Figure  12  (c)  it  is  apparent  that  a  straight  line  will  not  describe  the  curve 
accurately.  The  Hollomon  power  equation  has  the  form: 

a  =  K  em  or:  log  a  =  logK  +  m  log  € 

where  K  is  the  strength  coefficient  and  m  is  the  strain  hardening  exponent.  By  simply  doing 

a  least  squares  fit  of  log  a  as  a  function  of  log  e,  m  is  the  slope  of  the  line  and  the  intercept 
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is  logK. 

Using  this  method  for  the  solution  of  the  Hollomon  power  equation  coefficients  and  the 
method  previously  described  for  the  solution  of  the  Voce  coefficients,  for  all  of  the  tests,  the 
results  are  as  listed  in  Table  3. 

For  all  data  sets,  the  Hollomon  power  equation  fit  the  data  with  correlations  ranging 
from  .9492  to  .9836  and  the  Voce  equation  fit  with  correlations  of  from  .9937  to  .9993  for 
the  same  regions  of  the  curves.  From  the  point  of  correlation  alone  the  Voce  equation  would 
appear  to  fit  the  data  much  better.  Looking  at  a  typical  plot  of  the  true  stress/true  strain 
data  with  the  predicted  curves  overlayed  as  shown  in  Figure  13,  the  Voce  equation  does  fit 
the  data  much  better.  Neither  of  these  equations  fit  well  in  the  region  of  low  strains  so  for 
comparative  purposes  the  fifth  data  point  is  arbitrarily  chosen  as  the  first  point  compared 
and  the  last  point  used  in  the  curve  fit  is  the  last  considered.  Starting  from  this  premise, 
the  percentage  error  of  the  predicted  curve  is  calculated  by  the  equation: 

Percent  Error=  [ABS[  o  -  o  c]/  c  ]  *  100 

For  the  Hollomon  power  equation  the  maximum  single  point  error  is  8.30  percent  and  the 
maximum  average  error  for  any  single  curve  is  3.80  percent.  The  fit  appears  to  be  worse 
at  low  temperatures  and  gets  better  at  intermediate  temperatures  (-1 17°C).  The  best 
average  error  for  the  Hollomon  power  equation  is  1 .58  percent  at  -1 17°C. 


39 


TABLE  3 


VOCE  AND  HOLLOMON  EQUATION  COEFFICIENTS 

Temp 

sample 

Vooe 

Voce 

Voce 

Hollomon 

Hollomon 

(°C) 

number 

S~  (ksi) 

So  (ksi) 

A  (in/in) 

K  (ksi) 

m 

-189 

46 

491.0986 

1 75.2839 

1 .360575 

300.9788 

.1523583 

-176 

28 

399.506 

165.2733 

1 .095739 

285.465 

.1692908 

-170 

52 

389.2834 

169.0419 

1 .1 00783 

281.1485 

.1604666 

-154 

34 

356.417 

157.7226 

1 .077566 

262.1936 

.1601815 

-150 

47 

380.1898 

156.2891 

1.176447 

262.168 

.157335 

-140 

26 

385.2167 

154.5706 

1.309195 

261.0666 

.1674174 

-133 

25 

298.439 

148.9738 

.7496496 

249.049 

.1572877 

-132 

37 

309.1749 

144.9128 

.8607823 

245.0456 

.1610462 

-124 

24 

337.1254 

156.5421 

1.194708 

244.7035 

.1436531 

-117 

19 

257.3777 

143.9943 

.5703971 

232.392 

.143561 

-110 

49 

282.4977 

145.829 

.8121898 

232.2654 

.1440436 

-101 

23 

309.5079 

145.3076 

.9910815 

236.4248 

.1517314 

-91 

22 

253.6235 

140.7949 

.655102 

223.21 69 

.1433506 

-80 

50 

265.0807 

138.1622 

.781093 

220.8118 

.1458434 

-71 

20 

308.4355 

140.5525 

1 .072554 

230.6979 

.1585906 

-56 

31 

293.8192 

137.9099 

1 .065468 

222.25 

.1528224 

-30 

51 

279.0695 

132.6547 

1.01008 

214.241 

.15167 

0 

30 

271 .9845 

133.1512 

.9848082 

211.7802 

.1474638 

18 

36 

258.9012 

131.1185 

.9638342 

203.3129 

.1383482 

18 

35 

272.7748 

129.2532 

1.144794 

202.0873 

.1440793 

37 

27 

245.5962 

125.0026 

.8919725 

197.183 

.1418691 

The  Voce  equation  on  the  other  hand  has  a  worst  single  point  error  for  all  of  the  data 
sets  of  4.55  percent  and  a  worst  average  error  of  1 .21  percent.  The  best  average  error  was 
0.43  percent  at  -150°C.  The  maximum  errors  and  average  error  are  less  for  the  Voce 
equation  for  all  data  sets.  Sample  number  50,  shown  in  Figure  13,  has  the  smallest 
difference  between  average  errors  of  the  Hollomon  power  equation  and  Voce  equation  for  the 
tests  conducted.  Based  on  these  results  the  Voce  equation  was  determined  to  fit  the  data 
better,  and  the  Hollomon  equation  was  not  considered  for  further  investigation. 
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Figure  13 
Comparison  of  the  Voce  and  Power  Equations 
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2.    The  Voce  Equation  as  a  Function  of  Temperature 

Several  forms  of  equations  were  considered  for  a  fit  of  the  Voce  stress  coefficients  as 
a  function  of  temperature  including  exponential  and  logarithmic  forms.  The  final  form 
decided  upon  was  an  exponentially  decaying  function  typically  used  to  describe  a  thermally 
activated  event  such  as  overcoming  the  Peierls-Nabarro  stresses.  The  Voce  stress 
coefficients  as  a  function  of  temperature  are  found  to  be: 

Soo(T)=  Co  +  [Crexp(Tc/T(°K))] 
So(T)=  Co  +  [C1*exp(Tc/T(°K))] 
where  Co  and  C1  are  constants  with  units  of  stress,  Tc  is  a  constant  with  units  of 
temperature  in  degrees  Kelvin,  and  T(°K)  is  the  temperature  of  the  test  in  degrees  Kelvin. 
As  the  test  temperature  becomes  very  large  compared  to  Tc,  the  exponential  becomes 
approximately  equal  to  1  and  the  sum  of  Co  and  C1  then  becomes  the  value  of  stress  which  is 
asymptotically  approached  at  high  temperatures. 

The  method  for  solving  this  equation  was  similar  to  that  of  the  Voce  equation  except  in 
this  case  the  stress  coefficient  is  the  dependent  variable  and  the  quantity  exp(Tc/T)  is  the 
independent  term.  The  value  of  Tc  was  incremented  from  1  to  300°K  and  the  most  linear 
value  of  Tc  was  the  one  from  which  the  slope,  C1 ,  and  the  intercept,  Co,  were  calculated.  The 
linearity  was  again  measured  by  the  correlation  coefficient.  The  plot  of  the  stress 
coefficient  data  points  and  the  curve  fit  to  them  is  shown  in  Figure  14. 

The  characteristic  strain,  A,  plotted  against  temperature  yielded  a  random  scatter  of 
data  points  with  no  apparent  correlation  as  shown  in  Figure  15.  The  average  for  all  of  the  A 
values  was  found  to  be  0.99375  and  the  median  was  1 .01008  so  the  mean  was  arbitrarily 
chosen  to  be  the  value  of  A  across  all  temperatures  (A(T)=  constant=  .99375). 
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HSLA-100  True  Stress/True  Strain  -  Soo  vs.  Temperature 
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Figure  14 
Voce  Stress  Coefficients  as  a  Function  of  Temperature  (1) 
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HSLA-100  True  Stress/True  Strain  -  A  vs.  Temperature 
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Figure  15 
Characteristic  Strain  as  a  Function  of  Temperature 

By  combining  the  individual  equations,  the  Voce  equation  as  a  function  of  temperature 

and  strain  was  derived. 

S=  Soo(T)  -  {[Soo(T)  -  So(T)]*exp(-€/A)} 
where  S°o(T)=  1 91 .1 32  +  [39.1 0577*exp(1 70/T(°K))] 

So(T)=  -5397.923  +  [5509.1 05*exp(1/T(°K))] 
and  A=  .9937533 

A  test  expected  to  have  a  good  fit  and  another  test  expected  to  have  a  poor  fit  based  on  the 
differences  between  the  actual  S°°  and  So  and  that  predicted  by  the  above  equations  were 
selected  to  examine  the  accuracy  of  the  overall  predictions  based  on  these  temperature 
dependent  equations.  The  fit  of  the  predicted  curve  for  sample  46,  which  had  actual 
coefficients  right  on  the  predicted  curves,  and  sample  number  26,  which  has  S°°  predicted 
about  50  ksi  lower  than  the  test  indicated,  are  shown  in  Figure  16.  For  orientation  the 

label,  Voce(T),  is  the  curve  fit  as  a  function  of  e  and  T,  the  other  is  the  Voce  equation  curve 
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fit  for  that  individual  test.  The  characteristic  strain  and  the  S«>  have  some  sort  of  balancing 
relationship.  That  is,  even  though  a  difference  may  exist  between  the  predicted  values  of  S°o 
and  A  with  respect  to  the  value  of  the  constants  determined  from  the  individual  test,  the 
overall  stress-strain  prediction  may  be  excellent.  This  will  not  always  be  the  case, 
however. 

The  next  step  in  the  analysis  was  to  assume  that  the  characteristic  strain,  A,  is  a 
constant  as  already  assumed  and  recalculate  the  stress  coefficients.  Since  the  mean  and 
median  bracket  the  value  of  1 .0,  A  is  assumed  constant  and  equal  to  1 .0.  With  this 
assumption  the  values  for  S°°  and  So  are  recalculated  with  a  simple  slope  and  intercept 
calculation;  stress  the  dependent  variable,  strain  the  independent  variable,  S°°  equal  to  the 
intercept  and  the  slope  equal  to  -(S°°-So).  The  correlation  coefficient  for  all  version  2 
individual  curve  fits  were  greater  than  .990.  The  results  were  then  plotted  and  an 
exponential  curve  fit  to  to  the  new  stress  coefficient  data  points.  The  results  are  shown  in 
Figure  17;  by  the  correlation  coefficient  alone  these  curves  fits  the  new  coefficient  data 
better  than  those  initially.  The  coefficients  for  So  did  not  change  significantly  but  the  S°° 
coefficients  are  much  closer  to  the  predicted  curve,  increasing  the  correlation  coefficient  by 
greater  than  0.06.  In  the  same  manner  as  before,  the  Voce  coefficients  as  a  function  of 
temperature  are  substituted  into  the  Voce  equation  forming  an  equation  for  stress  as  a 
function  of  strain  and  temperature: 

S=  Soo(T)  -  {[Soo(T)  -  So(T)]*exp(-  e  /A)} 
where  S~(T)=  -78.8051  +  [292.1 489*exp(46/T(°K))] 

So(T)=  -5148.559  +  [5261.393*exp(1/T(°K))] 
and  A=  1 .0 

When  T(°K)  is  much  greater  than  Tc,  S~  approaches  213  ksi  and  So  approaches  113  ksi. 
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Figure  16 
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HSLA-100  Voce  Coefficient  S°o  vs  Temperature/A=1 
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Figure  17 
Voce  Stress  Coefficients  as  a  Function  of  Temperature  (2) 
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Fitting  the  new  Voce  equation  to  the  same  two  curves  as  before  gives  apparently  more 
satisfactory  results.  A  comparison  of  the  errors  of  the  two  versions  is  presented  in  Table  4. 
The  initial  coefficients  calculated  are  sometimes  better  than  those  calculated  for  version  2, 
but  overall  version  2  is  much  better.  Version  2  shows  generally  lower  single  point 

maximum  errors  and  lower  averages.  Figure  18  shows  the  modified  version  of  a(  e,  T) 
using  the  new  coefficients  as  a  function  of  temperature  overlayed  with  the  individual  test 
predictions  and  data  for  several  tests.  The  same  tests,  numbers  26  and  46,  are  presented  as 
are  two  other  tests,  samples  51  and  20,  the  tests  with  the  smallest  and  largest  average 
errors,  0.560  and  2.76  percent  respectively. 

TABLE  4 
COMPARISON  OF  ACCURACY  FOR  VOCE(e.T)  VERSIONS 


sample 

Temp 

Max  Error 

Avg  Error 

Max  Error 

Avg  Error 

number 

f°C) 

Version  1 

Version  1 

Version  2 

Version  2 

46 

-189 

1 1 .79377 

8.095257 

2.179396 

.8623438 

28 

-176 

7.796861 

5.141882 

3.878173 

1 .28702 

52 

-170 

3.967472 

2.03051 

2.997622 

.9321079 

34 

-154 

3.87577 

1 .433683 

3.495088 

1.046242 

47 

-150 

2.674214 

1 .028046 

3.1547 

1 .370743 

26 

-140 

5.042438 

2.307615 

4.677677 

2.134145 

25 

-133 

3.60902 

1.144696 

3.60868 

.9944419 

37 

-132 

6.272017 

.7508031 

6.287417 

1.196656 

24 

-124 

5.571671 

2.8092 

5.203831 

2.313026 

19 

-117 

4.789085 

1 .638871 

6.210592 

2.020577 

49 

-110 

3.813052 

1.196474 

4.02332 

1 .452825 

23 

-101 

4.538417 

3.012768 

3.873194 

2.032468 

22 

-91 

3.561259 

1 .429343 

5.00009 

1 .782425 

50 

-80 

5.269707 

1.066175 

5.666542 

1 .546789 

20 

-71 

5.652589 

3.804017 

4.19193 

2.764634 

31 

-56 

3.811433 

2.316394 

2.7988 

1 .337599 

51 

-30 

2.569845 

1 .258869 

2.944652 

.5596738 

30 

0 

3.998918 

2.392466 

3.413519 

1 .762665 

36 

18 

2.750092 

.8317102 

2.195829 

.7063553 

35 

18 

2.871258 

.8925604 

3.431084 

1 .259844 

27 

37 

5.589019 

1.994318 

6.239485 

2.401873 
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Figure  18 
Final  Voce(  e,T)  Comparison  Plots 
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Looking  at  the  Voce  equation  fit  to  the  data,  there  appears  to  be  three  regions  where 
the  fit  is  poor.  The  low  strain  region  of  all  tests  is  generally  the  region  where  the  fit  is  the 
worst;  the  Voce  curve  always  overestimating.  The  second  region  where  the  curve  fit  is  poor 
is  at  the  maximum  load  point  where  the  stress  is  generally  predicted  lower  than  the 
observed  values,  the  last  region  is  the  very  high  strain  region  of  the  curve  near  failure.  In 
this  region  the  prediction  is  higher  than  actual  because  of  the  rapid  drop  in  load  as  the 
specimen  fails.  It  is  not  possible  to  just  change  the  value  of  S°°  or  So  to  an  observed  value 
without  changing  the  value  of  the  characteristic  strain  because  of  the  close  interaction  of  the 
coefficients.  It  is  possible  to  obtain  better  results  in  all  of  these  areas  by  splitting  the 
stress/strain  curve  into  two  regions,  one  from  the  yield  point  to  the  maximum  load  point 
and  another  from  the  maximum  load  point  to  failure.  Because  of  the  desire  to  obtain  a  single 
equation  describing  stress  as  a  function  of  strain  and  temperature  this  path  was  not  pursued 
in  great  detail.  However  the  data  suggests  that  there  may  be  a  change  in  the  controlling 
strain  hardening  mechanism  near  the  point  of  maximum  load. 

D.   BRIDGMAN  CORRECTED  TRUE  STRESS  VERSUS  TRUE  PLASTIC  STRAIN 

The  Bridgman  correction  is  a  stress  correction  factor  which  corrects  for  the  triaxial 
stress  state  that  is  generated  in  the  neck  of  a  tensile  sample.  Since  the  necking  phenomenon 
begins  at  approximately  the  point  of  maximum  load  the  correction  becomes  greatest  in 
regions  past  the  maximum  load  point.  In  the  derivation  of  the  correction  factor,  Bridgman 
was  able  to  describe  the  correction  in  terms  of  some  relatively  easy  to  measure  parameters 
of  the  tensile  specimen,  the  diameter  of  the  specimen  at  the  minimum  diameter  and  the 
radius  of  curvature  in  the  neck.  The  Bridgman  correction  factor  is  given  as: 
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a  /aAv=  1/[(1  +  AR^JD)  Ln(1  +  D/4RC)] 

where  a  is  the  Bridgman  corrected  stress  and  is  always  lower  than  the  measured  stress,  o~ 
Av.  Rc  is  the  radius  of  curvature  of  the  neck  and  D  is  the  minimum  diameter  of  the  neck. 

[Ref.  13] 

Sample  number  43  was  interrupted  during  testing  at  several  intervals  past  the 
maximum  load  point  and  the  instantaneous  diameter  and  radius  of  curvature  measured. 
Several  methods  for  measuring  the  radius  of  curvature  were  tried  including  projecting  the 
sample  image  and  using  a  circle  template  to  find  a  radius  proportional  to  the  known  diameter 
of  the  specimen.  Another  method,  the  use  of  a  machinists  radius  gage  set  was  less  successful 
because  the  gage  could  not  fit  between  the  grips.  The  most  successful  method,  especially  at 
low  temperatures  was  to  use  a  drill  bit  index.  By  back  lighting  the  specimen  the  best  radius 
fit  could  be  found  at  temperatures  above  0°C.  At  low  temperatures  the  frost  which  formed 
on  the  specimen  melted  when  touched  by  the  drill  bit,  so  the  largest  bit  which  melted  all  the 
frost  through  the  neck  was  the  diameter  of  the  neck.  Drill  bits  incremented  in  1/64  inch 
gives  radii  accurate  to  1/128  of  an  inch  or  about  .008  inches.  Another  advantage  of  the 
drill  index  method  is  that  it  is  not  as  subjective  as  the  projection  method,  the  drill  index  can 
measure  with  a  smaller  increment  size  than  a  radius  gage  set  again  making  it  more  accurate. 

The  plot  of  all  of  the  Bridgman  correction  data  collected  and  the  equation  of  an  computer 
generated  polynomial  fit  are  shown  in  Figure  19.  Since  the  tests  were  conducted  at  a  constant 
displacement  rate,  the  diameter  change  is  constant  and  the  change  in  radius  of  curvature 
must  be  responsible  for  the  more  rapid  changes  in  the  Bridgman  correction  near  the 
maximum  load  point  and  again  near  failure. 

Applying  the  correction  to  the  actual  data  was  accomplished  by  knowing  the  maximum 
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load  for  a  particular  test  and  the  strain  at  that  particular  point;  all  strains  past  the  necking 
strain  were  calculated  and  the  Bridgman  correction  calculated  from  the  equation  shown  in 
Figure  1 9.  For  strains  prior  to  necking  the  diameter  was  changing  and  the  radius  of 
curvature  was  assumed  to  remain  constant.  Although  this  is  known  not  true  in  fact,  the 
correction  factor  remains  relatively  constant  for  increases  in  radius  of  curvature  so  this  is 
a  good  approximation  prior  to  necking.  Converting  the  true  strain  into  true  plastic  strain  is 
accomplished  by  subtracting  the  elastic  strain  from  the  total  strain. 
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Figure  19 
Collection  of  all  Bridgman  Correction  Data  Points  for  HSLA-100  Tests 


Applying  these  corrections  to  an  actual  test  yielded  results  shown  in  Figure  20.  On  these 
plots  the  change  in  slope  at  the  maximum  load  point  is  much  more  obvious.  On  the 
logarithmic  plot,  there  appears  to  be  two  separate  and  overlaying  exponential  plots 
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intersecting  at  the  maximum  load  point.  Toward  the  high  strain  regions  of  the  curve  the 
stress  is  dropping  off,  showing  as  a  hook  in  the  logarithmic  plot,  another  unexpected  result. 
For  the  single  Voce  equation  fit  for  this  curve,  the  correlation  was  .97;  a  large  decrease  in 
the  accuracy  of  the  prediction. 

The  Bridgman  correction  factor  for  all  tests  are  close  to  the  prediction  curve  used  so  the 
confidence  in  the  Bridgman  Correction  factor,  as  it  was  predicted  and  used,  is  high.  The  true 
plastic  strain  is  too  straightforward  to  account  for  the  deviations  observed.  There  was  no 
observed  reason  during  the  testing  which  could  describe  the  rapid  change  in  stress  carrying 
ability  at  the  maximum  load  point  and  then  again  near  failure.  The  Bridgman  correction  was 
not  given  any  further  consideration  because  it  is  felt  that  it  is  not  properly  correcting  for 
the  triaxiality  effects  on  the  true  stress  of  this  material. 

E.    MECHANICAL  PROPERTIES  OF  HSLA-1 00  WITH  TEMPERATURE 

Ductility,  0.2  percent  yield  strength,  and  ultimate  tensile  strength  are  more 
conventional  material  properties  specified  in  metals  specifications.  For  this  HSLA-1 00  the 
ultimate  tensile  strength  was  to  be  recorded  for  information  purposes  only,  the  0.2  percent 
offset  yield  strength  was  to  be  between  100  and  120  ksi,  and  the  minimum  reduction  in  area 
was  to  be  45  percent.  The  HSLA-1 00  used  in  this  research  showed  many  of  the  expected 
trends  with  regards  to  temperature.  A  tabular  compilation  of  these  results  are  provided  in 
Table  5. 

1.    0.2  Percent  Offset  Yield  Strength 

The  0.2  percent  offset  yield  strength  is  one  characteristic  which  showed  a  wide 
degree  of  variability  in  testing.  The  results  of  all  tests  are  shown  in  Figure  21 .  In  most 
tests,  the  yielding  behavior  of  this  steel  was  such  that  load/displacement  curve  rounded  off 
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slightly  as  the  yield  point  was  approached  and  then  leveled  off  to  horizontal  for  a  short 

period  in  what  may  have  been  a  second  yield  point  for  the  material. 
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TABLE  5 
MATERIAL  PROPERTIES  OF  HSLA-100 


Ductility 
AA/A 

.5215 

.5977 

.5701 

.6217 

.599 

.6551 

.6168 

.6198 

.6799 

.6578 

.6467 

.6643 

.6643 

.6794 

.6670 

.6766 

.7173 

.7095 

.6898 

.7130 

.7037 


The  ductility,  measured  by  reduction  in  area  (AA/A),  of  the  HSLA-100  is  shown  in 
Figure  22.  The  ductility  was  lowest  at  low  temperatures  with  a  minimum  reduction  in  area 
of  52  percent. 

3.    Ultimate  Tensile  Strength 

The  true  stress  at  the  maximum  load  point  is  shown  in  Figure  23  as  a  function  of 
temperature.    This  data  exhibit  the  expected  changes  with  temperature,  increasing  tensile 
strengths  as  the  temperature  decreases. 
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Figure  21 
0.2  Percent  Offset  Yield  Strength 
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HSLA-100  Ductility  vs.  Temperature 
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Figure  22 
Ductility  of  HSLA-100 
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HSLA-100  Ultimate  Tensile  Strength 
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Figure  23 
HSLA-100  Ultimate  Tensile  Strength 


F.   FRACTURE  CHARACTERISTICS  AND  METALLOGRAPHY 

Fractographic  studies  of  the  macroscopic  and  microscopic  features  of  the  fracture 
surfaces  were  conducted  using  light  microscopy  and  scanning  electron  microscopy. 

The  most  notable  macroscopic  feature  of  the  this  material  was  the  sequence  of  events 
during  the  failure  process.  Stated  in  a  straightforward  manner,  only  two  of  the  1 1  samples 
tested  at  temperatures  from  -56°C  to  -140°C  fractured  into  two  pieces  during  actual 
testing.  In  this  temperature  range  and  below  the  first  fracture  was  a  vertical  split  or 
delamination  in  the  neck  and  running  for  a  total  length  of  about  1/4  of  an  inch.  A  photograph 
of  one  of  these  splits  is  shown  in  Figure  24.  In  a  few  of  the  samples  there  were 
perpendicular  vertical  splits,  but  in  all  cases  there  was  one  major  split.  For  all  of  these 
samples  there  was  the  loud  pop  associated  with  a  sample  failure  but  the  sample  remained 
intact.  The  deformation  in  this  temperature  range,  at  this  point,  is  presumed  to  have 
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occurred  in  a  jump  since  the  tests  were  terminated  by  the  MTS  controller  due  to  the  error 
between  control  and  actual  displacement  voltages  becoming  too  large.  Metal  log  raphic 
examination  of  the  split  shows  that  it  is  parallel  to  the  banding  noted  in  the  microstructure. 
Figure  25  is  a  photograph  showing  the  microstructure  and  its  orientation  with  the  split. 

From  the  geometry  of  the  sample,  all  plastic  deformation  is  restricted  to  a  small  length 
in  the  center  of  the  sample.  The  actual  length  of  localized  deformation  depends  on  the  yield 
strength,  strain  hardening,  ultimate  tensile  strength  and  initial  diameter  of  the  sample. 
With  the  radius  of  curvature  of  the  hourglass  portion  equal  to  1 .1  inches,  the  plastic 
deformation  occurs  in  the  center  0.20  to  0.25  inches  of  the  initial  sample.  Looking  at  the 
sample  and  measuring  from  the  base  of  the  hourglass  portion  (the  elastic  region  of  the 
sample)  the  splits  appear  to  end  at  about  the  location  of  the  boundary  between  elastically  and 
plastically  strained  material. 

Fracture  profiles  of  Figure  26  show  some  specimens  tested  in  the  transition  region  of 
-140°C  to  -189°C.  The  fracture  surfaces  perpendicular  to  the  load  of  all  tests  except  that 
at  -189°C  appear  to  have  occurred  at  about  the  location  of  the  elastic-plastic  boundary.  At 
this  point  it  is  also  noted  for  samples  26  and  47  in  Figure  26  (b)  (-140°C  and  -150°C) 
that  there  was  failure  both  in  and  outside  of  the  neck. 
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Figure  24 

HSU\-100  Sample  19  (-117°C)  Vertical  Split 

1  Ruler  Division  =  1/64  inch 


Figure  25 
HSLA-100  Sample  52  (-170°C)  Cross  Section  of  Neck 
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(a.) 


(b.) 

Figure  26 
Fracture  Profiles 
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Scanning  electron  microscopy  observations  on  the  fracture  surfaces  shows  many 
important  features  which  can  be  linked  to  the  fracture  profiles.  The  various  viewsof  the 
fracture  surfaces  are  shown  in  Figure  27.  At  temperatures  from  room  temperature  to 
-30°C  the  typical  cup  and  cone  of  ductile  fracture  were  noted,  with  the  appearance  of  early 
stages  of  delamination.  The  higher  magnification  photographs  show  some  of  the  sites  where 
microvoid  coalescence  was  initiated.  At  -71  °C  the  area  of  delamination  has  opened  up  and 
the  cup  and  cone  surface  is  no  longer  evident.  At  -1 01  °C  the  surface  has  taken  on  a  jagged 
appearance  and  the  high  magnification  photograph  shows  a  still  predominantly  ductile 
fracture  surface  with  small  regions  of  quasi-cleavage.  At  -140°C  and  -150°C  the  mode  of 
failure  was  dependent  on  the  the  location  of  the  fracture  surface  relative  to  the  neck.  In  the 
areas  where  the  failure  was  out  of  the  neck,  the  fracture  surface  was  brittle  (cleavage)  and 
in  the  areas  where  it  was  in  the  neck  ,  the  surface  was  all  ductile.  Below  this  temperature 
all  of  the  fracture  surfaces  were  out  of  the  neck,  except  for  the  test  at  -1 89°C.  At  -1 89°C 
there  also  appeared  to  be  more  ductile  areas  on  the  fracture  surface  than  the  surfaces  at 
immediately  higher  temperatures  although  the  fracture  was  still  predominantly  brittle. 
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Figure  27(a.) 
Sample  30  (0°C)  Fracture  Surface,  Low  Magnification 
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Figure  27(b.) 
Sample  51  (-30°C)  Fracture  Surface,  Low  Magnification 
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Figure  27(c) 
Sample  30  (0°C)  Fracture  Surface,  High  Magnification 


Figure  27(d) 
Sample  20  (-71  °C)  Fracture  Surface,  Low  Magnification 
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Figure  27(e.) 
Sample  23  (-101°C)  Fracture  Surface,  Low  Magnification 
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Figure  27(f.) 
Sample  23  (-101°C)  Fracture  Surface,  High  Magnification 
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Figure  27(g.) 
Sample  26  (-140°C)  Fracture  Surface,  Low  Magnification 
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Figure  27(h.) 
Sample  47  (-150°C)  Fracture  Surface,  Low  Magnification 
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Figure  27(i.) 
Sample  52  (-170°C)  Fracture  Surface,  High  Magnification 


67 


V.  DISCUSSION 

A.   MODIFICATIONS  /  IMPROVEMENTS  IN  THE  USE  OF  THE  VOCE  EQUATION 

In  the  derivation  of  the  constituitive  Voce  equation,  a(€,T),  the  characteristic  strain,  A, 
was  taken  to  be  a  constant  equal  to  1 .0.  This  value  is  the  average  of  all  of  the  values  of  A 
obtained.  Some  of  the  variables  of  testing  which  affect  its  value  are:  which  region  of  the 
stress-strain  curve  is  being  evaluated;  how  far  the  test  was  run  toward  actual  failure;  and 
how  large  a  deviation  existed  in  the  curve  about  the  maximum  load  point.  There  is  also  an 
intimate  relationship  between  S°°,  So  and  A  so  that  if  one  is  changed  the  others  also  have  to 
be  changed  for  the  curve  to  fit  the  data.  The  characteristic  strain  is  believed  to  be  a  constant 
for  this  series  of  tests  but  may  be  dependent  on  the  strain  rate.  For  this  series  of  tests  the 
strain  rate  was  relatively  low,  about  1 .0x1 0"^  in/in-sec.  At  higher  strain  rates  the  value 
of  S«>  would  be  expected  to  raise,  increasing  the  quantity  Ln(S°°  -  S)  for  a  given  strain. 
The  overall  effect  would  be  to  lower  the  slope  of  the  strain-Ln(S°°  -  S)  curve,  which  is  the 
characteristic  strain. 

There  are  two  definite  strain  hardening  regions  for  this  material  with  a  transition 
occurring  at  about  the  maximum  load  point.  In  an  analysis  of  the  two  regions,  the  average 
value  of  A  from  the  yield  point  to  necking  was  0.1829  and  the  average  from  the  maximum 
load  point  to  failure  was  1 .522  showing  a  large  difference  between  the  two  regions.  It  is  for 
this  reason  that  the  characteristic  strain  should  not  be  dropped  from  the  Voce  equation 
describing  the  stress  as  a  function  of  strain  and  temperature.  Further,  this  suggests  that 
future  research  should  attempt  to  determine  if  there  is  a  change  in  the  strain  hardening 
mechanism  after  necking  occurs. 
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Using  separate  versions  of  the  Voce  equation  to  describe  the  two  strain  hardening  regions 
seems  initially  to  more  accurately  fit  the  curves  with  correlations  of  .999  occurring 
frequently.  In  order  to  effectively  use  a  two  region  /  two  equation  form  of  solution,  the 
transition  point  must  be  accurately  predicted  also.  For  the  case  of  the  Voce  equation  and 
HSLA-100  this  point  is  the  maximum  load  point.  Reference  12  points  out  that  the  slope  of 
the  stress/strain  curve  at  the  maximum  load  point  equals  the  true  stress  at  that  point.  Based 
on  this  fact  and  using  the  derivative  of  the  Voce  equation  for  the  conditions  of  maximum  load 
the  stress  at  maximum  load  is  described  by  the  equation: 

da/de  =  a  max  )oad  =  S~/(1  +  A) 

and  for  the  case  of  A=  1  the  stress  at  maximum  load  =  S<*V2.  The  stresses  calculated  using 
this  equation  and  the  actual  ultimate  tensile  strengths  found  experimentally  are  plotted  in 
Figure  28.  The  individual  Voce  equation  coefficients  for  the  test  and  the  general  Voce 
coefficients  as  a  function  of  temperature  are  all  significantly  lower  than  the  actual  values 
with  the  worst  results  near  room  temperature. 

From  the  experimental  results  it  is  known  that  the  predicted  value  of  stress  at  the 
maximum  load  point  is  always  low  but  this  also  raises  a  question  of  confidence  in  the  Voce 
equation  for  describing  actual  material  behavior.  The  Voce  equation  does  fairly  accurately 
describe  the  stress/strain  behavior  of  HSLA-100  but  there  are  material  inconsistencies 
observed  and  not  fully  explained  by  this  equation.  If  this  problem  of  the 
transition/maximum  load  point  can  be  resolved,  then  the  Voce  equation  or  one  of  its 
derivatives  could  be  used  very  well  to  predict  the  behavior  of  HSLA-1 00  in  terms  of  true 
stress/true  strain.  For  one  test,  the  two  region  Bridgman  corrected  curve  was  fit  and  is 
shown  in  Figure  29.  The  results  appear  to  be  very  good  except  at  the  very  end  in  the 
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vicinity  of  the  hook,  where  the  sample  is  starting  to  fail. 

B.   SAMPLE  GEOMETRY  RELATIONSHIP  TO  FRACTURE 

The  hourglass  specimen  was  used  to  localize  the  deformation  in  the  specimen  enabling  the 
true  stress  state  to  be  accurately  determined  through  necking  up  to  fracture.  Fracture 
characteristics  such  as  splitting  only  in  the  plastic  region  and  final  fractures  near  the 
elastic-plastic  interface  bring  rise  to  questions  of  the  applicability  of  information  obtained 
from  hourglass  specimen  tests.  However,  preliminary  tensile  tests  using  uniform  gage 
samples  at  David  Taylor  NSRDC  have  fractured  in  an  identical  fashion  as  these  hourglass 
specimens. 


HSLA-100  Prediction  of  True  Stress  at  Maximum  Load 
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Figure  28 
Prediction  of  the  Ultimate  Tensile  Strength 
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Figure  29 
Two  Curve  Fit  for  the  Bridgman  Corrected  True  Stress  /  True  Plastic  Strain 
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The  aspect  of  the  fracture  that  is  most  puzzling  is  the  fact  that  all  cleavage  fractures 
occurred  outside  of  the  neck  where  the  nominal  stress  is  much  lower  than  in  the  neck.  It 
appears  that  the  large  plastic  strain  experienced  within  the  necked  region  has  raised  the 
cleavage  stress.  All  fractures  in  the  neck  occurred  by  microvoid  coalescence  while  all 
fractures  outside  the  neck  occurred  by  cleavage.  Some  samples  in  the  transition 
temperature  range  exhibited  regions  of  ductile  fracture  in  the  neck  and  brittle  fracture 
outside  the  neck.  The  effects  of  large  prestrains  on  cleavage  fracture  should  be  the  topic  of 
future  research. 

The  rapid  increase  in  strength  with  decreasing  temperature  indicates  that  the 
Peierls-Nabarro  stress  is  an  important  component  of  the  strength  of  this  alloy,  especially 
at  low  temperatures.  Thus,  the  exponential  equation  describing  the  temperature  dependence 
of  the  strength  coefficients,  So  and  S°°,  in  the  Voce  equation  should  be,  at  least,  correct  in 
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form.  The  characteristic  strain,  A,  in  the  Voce  equation  is  a  measure  of  the  strain  hardening 
rate  of  the  material.  Considerable  scatter  in  the  value  for  A  was  found,  but  no  evidence  of 
any  temperature  dependence  of  A  was  observed.  This  indicates  that  the  mechanism 
controlling  the  strain  hardening  in  this  alloy  is  probably  independent  of  temperature. 
However,  the  fact  the  stress-strain  curve  has  what  appears  to  be  different  strain  hardening 
regions  before  and  after  the  maximum  load  may  indicate  that  the  controlling  mechanism 
changes  once  a  critical  strain  is  exceeded.  If  this  is  the  case  then  a  single  constituitive 
equation  like  the  Voce  equation  will  need  to  be  modified  to  reflect  the  change  in  the  strain 
hardening  mechanism.  However,  the  single  Voce  equation  which  incorporates  the  effect  of 
temperature  appears  to  produce  a  very  adequate  description  of  the  stress-strain  behavior  of 
this  alloy  at  this  quasi-static  strain  rate. 
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VI.    CONCLUSIONS  AND  RECOMMENDATIONS 

The  Voce  equation  is  quite  accurate  in  describing  the  the  tensile  behavior  of  HSLA-100 
steel  but  has  some  areas  of  weakness.  In  low  strain  regions  (less  than  about  0.01 )  the  Voce 
equation  always  overestimates  the  actual  stress.  The  Voce  equation  also  underestimates  the 
stress  at  the  point  of  maximum  load. 

Taken  as  a  function  of  temperature,  the  value  of  A  appears  to  be  a  constant;  however  the 
temperature  dependence  of  So  and  S°o  are  accurately  described  by  an  exponential  function  of 
the  general  form: 

S=Co  +  [C1*exp(Tc/T)] 
where  Co,  C1 ,  and  Tc  are  constants  and  T  is  the  temperature  of  the  test  in  degrees  Kelvin 
(°K).  Using  least  squares  curve  fit  of  observed  data  to  the  equation,  the  coefficients  are 
defined  as  follows: 

So=  -5148.559  +  [5261.393  *  exp(1/T)]  (ksi) 

Soo=  -78.8051 1  +  [292.1489  *  exp(46/T)]  (ksi) 

A=  constant=  1 .0  (in/in) 
These  equations  can  be  substituted  into  the  Voce  equation  to  form  a  constitutive  equation  for 
the  true  stress  in  HSLA-100  at  low  strain  rates  as  a  function  of  strain  and  temperature: 

a(  €,T)=  Soo(T)  -  [(Soo(T)  -  So(T))exp-(e/A)]. 
Recommendations  for  further  study  include  evaluating  the  effect  of  higher  strain  rates 
on  the  stress  coefficients  of  the  Voce  equation  with  particular  attention  to  the  characteristic 
strain.  A  data  acquisition  system  for  such  experiments  should  in  addition  to  monitoring  the 
load,  should  simultaneously  be  able  to  measure  the  the  detailed  plastic  flow  in  the  region  of 
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the  neck.  The  use  of  a  bit  mapping  scheme  or  diffraction  of  laser  beams  may  work  well.  This 
is  important  to  be  able  to  more  accurately  predict  the  actual  triaxial  stress  correction  and 
the  corresponding  stress  state  in  the  neck.  The  mechanism  of  deformation  and  the  effect  of 
prestrain  on  cleavage  fracture  also  need  further  study. 

Another  topic  poorly  understood  is  the  effect  of  the  delamination  or  splitting,  on  the 
stress  state  in  the  necked  region  of  the  test  specimen.  Cleavage  fractures  always  initiated 
transverse  to  the  loading  direction  at  the  end  of  the  axial  split,  a  region  of  much  lower  stress 
than  in  the  neck  for  these  hourglass  samples.  Why  does  the  material  exhibit  cleavage  here 
instead  of  in  the  neck  where  the  stresses  are  much  larger  ? 
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APPENDIX  A 

INTERIM  SPECIFICATION  FOR  TRIAL  COMMERCIAL  PRODUCTION  OF  HSLA-100  STEEL 

PLATES 

Melting.  Refining  and  Casting 

The  heat  shall  be  fully  killed  and  produced  to  fine  grain  practice.  It  shall  be  made  with  a 
low  sulfur  practice,  vacuum  degassed  and  argon  injected  with  CaSi  or  Mg  for  sulfide  shape 
control.  The  heat  shall  be  ingot  cast  with  bottom-pour  molds  to  insure  good  surface. 
Chemical  Composition 


The  Chemical  composition  shall  be  as  shown  in  Table  I. 

Taraet  for 

Max.  %  bv  Weiaht  Unless 

Element 

First  Heat 

a  Range  is  Indicated 

Carbon 

0.04 

0.06 

Manganese 

0.90 

0.75  - 1 .05 

Phosphorus 

ALAP* 

0.015 

Sulfur 

ALAP* 

0.006 

Silicon 

0.25 

0.40 

Nickel 

3.5 

3.35  -  3.65 

Chromium 

0.60 

0.45  -  0.75 

Molybdenum 

0.60 

0.55  -  0.65 

Copper 

1.60 

1 .45  -  1 .75 

Columbium 

0.025 

1 .45  -  1 .75 

Aluminum 

0.030 

0.020  -  0.040 

Nitrogen 

0.010 

0.015 

*  As  Low  As  Possible 
Hot  Rolling 

Plates  1/4, 3/4, 1-1/4,  and  2  in.  thick  shall  be  rolled.  Extra  care  shall  be  taken  to 
minimize  rolled-in  scale  that  could  interfere  later  with  achieving  an  adequate  cooling  rate 
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during  quenching  from  the  solution  treating  temperature.  The  plates  shall  be  roller  leveled 
while  still  warm  after  rolling. 
Heat  Treatment 

All  of  the  plates  shall  be  solution  heat  treated  for  one  hour  at  1650°F  and  platen 
quenched  with  high  pressure  water  jets  from  above  and  beneath  the  plate.  The  quench  water 
shall  not  exceed  100°F  to  ensure  an  efficient  quench. 

The  plates  shall  be  given  an  age  hardening  treatment  using  temperatures  and  times 
determined  for  each  plate  by  preliminary  tensile  testing  of  samples  from  cupons  aged  at 
various  conditions.  Aging  conditions  for  the  plates  shall  be  chosen  so  as  to  give  the  tensile 
properties  listed  in  Table  II. 
Mechanical  Properties 

The  heat  treated  material  shall  meet  the  tensile  property  requirements  specified  in 
Table  II  and  the  impact  property  requirements  specified  in  Table  III. 

Table  II  Tensile  Properties 

Ultimate  Tensile  To  be  Recorded  for 

Strength,  psi    Information  Only 

Yield  Strength,  <0.75  in.  >Q.75  in. 

0.2%  Offset,  psi  1 00,000  to  1 00,000  to 

120,000  115,000 

Min.  Elongation  17  18 

in  2  in.,  % 

Min.  Reduction  in  -  -  45 

Area,  Round  Specimen,  % 

The  tensile  properties  shall  be  determined  as  the  average  value  of  duplicate  specimens 

from  each  plate  tested  in  accordance  with  ASTM  method  of  testing  E8.  Full  thickness  flat 

specimens  shall  be  tested  for  the  1/4-in.  thick  plate  and  standard  round  specimens  0.505 
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in.  diameter  shall  be  tested  for  the  plates  3/4  in.  thick  and  thicker.  All  specimens  shall  be 
taken  tranverse  to  the  primary  rolling  direction. 

Table  III  -  Impact  Properties 


Plate 

CVN 

Thickness 

Specimen 

Test 

Energy,3 

Test 

in. 

Size 
5  mm  x  10  mm 

Temp..  °F 
0±3 

ft -lb 

0.25 

28 

Charpy  V-Notch 

-120  ±3 

15 

Tranverse 

0.75, 1 .25, 

10  mm  x  10  mm 

0±3 

55 

2.00 

-120  ±3 

30 

a  Avg.  of  three  tests,  minimum. 

The  Charpy  impact  properties  shall  be  determined  in  accordance  with  ASTM  method  of 
testing  E23.  Three  tests  transverse  to  the  final  rolling  direction  of  the  plate  shall  be 
conducted.  No  single  value  shall  fall  below  the  minimum  average  specified  in  Table  III  by 
more  than  5  ft-lb  for  standard  specimens  and  2-1/2  ft-lb  for  half  size  specimens. 
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APPENDIX  B 
BASIC  Computer  Program  for  Calculation  of  Voce  Equation  Coefficients 


REM**  THIS  PROGRAM  IS  TO  TAKE  LOAD/DISPLACEMENT  DATA  FROM  A  DATA 

REM**  FILE,  PUT  IT  INTO  ARRAYS  AND  THEN  CALCULATE  THE  VOCE  EQUATION 

REM**  COEFFICIENTS.  VERSION  1.3  USES  TRUE  STRESS  VS  TRUE  STRAIN 

REM** 

REM**  This  program  is  written  for  the  Apple  Macintosh  PC  in  Microsoft  BASIC 

REM**  Version  3.0 

a*********************************************************** 

REFERENCE:  Fatigue,  Tensile,  and  Relaxation  Behavior  of  Stainless  Steels; 
by  J.B.  Conway,  R.H.  Stentz,  and  J.T.  Berling;  1975  by  Technical 
Information  Center,  Office  of  Information,  United  States  Atomic  Energy 

Commission. 

************************************************************ 

Major  variables  used  in  the  program: 

A1     Instantaneous  cross  sectional  area 

CORR     Local  maximum  correlation  for  a  pass 

CORRMAX    The  maximum  correlation  found  for  the  data  points 

DO     Initial  sample  diameter 

DISP    Diametral  displacement 

LOD     Load 

MAXLOD     Load  at  UTS 

MAXLODPT     Number  of  the  data  point  of  the  max  load 

R    Correlation 

SINF     Incremented  S-infinity 

SINFO     Value  of  S-infinity  for  the  maximum  correlation 

SO     Load  where  plastic  deformation  begins=  the  0.1  %  yield  strength 

STRAIN     True  strain 

STRESS     True  stress 

STRESSU     Ultimate  tensile  strength 

REM** 

REM**  S=  Soo  -(Soo  -  So)  eA-(strain/A) 

REM** 

REM**    THE  METHOD  OF  SOLVING  THIS  EQUATION  IS  AN  ITERATIVE  ESTIMATION 

REM**  OF  Soo  UNTIL  THE  PLOT  OF  Ln(Soo-S)  vs  STRAIN  (SEMILOG  PLOT)  IS  FOUND 

REM**  TO  BE  MOST  LINEAR  BY  THE  HIGHEST  CORRELATION  COEFFICIENT. 

REM**  AT  THIS  VALUE  OF  S°o,  A  IS  THE  SLOPE  OF  THE  LINE  AND 

REM**  THE  INTERCEPT  WILL  GIVE  THE  VALUE  OF  (Soo-So)  FROM  WHICH  So  CAN 

REM**  BE  CALCULATED 

REM 

REM**  DIMENSIONING  THE  CALCULATION  ARRAY  AND  INPUT  OF  DATA 
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DIM  LOD(500)  ,  DISP(500)  ,  STRESS(500)  ,  STRAIN(500) 

NEXTFILE: 

PRINT 

PRINT 

PRINT  "Input  the  name  of  your  input  file  name  with  the  load/displacement" 

PRINT  "data  in  it.  The  data  should  be  in  the  order  of  [load,  displacement  ] " 

PRINT 

PRINT  "INPUT  FILE  NAME:"; 

BEEP 

INPUT  DATS 
OPEN  DAT$  FOR  INPUT  AS  #1 
PRINT 
PRINT 
PRINT  "Input  the  sample  original  diameter  in  inches"; 

BEEP 

INPUT  DO 

R0=  1 .1         'Radius  of  curvature  of  the  hourglass  samples 

PRINT 

PRINT 
Pl=  3.14159 
N=0 

SMAX=0  :  MAXLOD=  0 

REM**    LOOP  THROUGH  TO  STORE  LOAD/DISP  DATA  IN  ARRAYS  FOR  FUTURE 
REM**     CALCULATIONS 

WHILE  NOT  EOF(1) 
N=N+1 

INPUT#1,LOD(N),DISP(N) 

IF  LOD(N)  >  MAXLOD  THEN     'Looking  for  the  point  where  necking  begins 
MAXLOD=  LOD(N) 
.  MAXLODPT=  N 
END  IF 

WEND 
CLOSE  #1 
REM 
AGAIN: 

PRINT  N;"  sets  of  data  have  been  read  from  file  ";  DAT$; "  stress/strain  and  " 
PRINT  "are  now  being  calculated." 
REM        Calculation  of  stresses  and  strains 

FOR  I  =  1  TO  N 
A1  =  (PI/4)*((D0  -  DISP(I))A2)         'instantaneous  cross  sectional  area 
STRAIN(I)=  2*LOG(D0/(D0-DISP(l)))  '2  x  diametrial  strain=  longitudinal  strain 
STRESS(I)=  LOD(l)/A1  'insantaneous  true  stress 

NEXT  I 

PRINT 

PRINT  "Necking  begins  at  point  ";MAXLODPT; "  for  load" 

PRINT 
BEGINNING: 
SMAX=  0 
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PRINT 

PRINT  "Enter  the  range  of  points  you  desire  to  look  at  using  data  point  numbers" 

PRINT "  separated  by  a  comma."; 

BEEP 
INPUT  STARTPT.STOPPT 
'**  Looking  for  the  max  stress  in  the  interval  we  are  looking  at 

FOR  1=  STARTPT  TO  STOPPT 

IF  STRESS(I)  >  SMAX  THEN  SMAX=  STRESS(I) 

NEXT  I 
■ 

PRINT 

PRINT  "Sums  of  the  data  are  now  being  calculated  for  a  best  correlation  fit." 

PRINT 

PRINT 

REM     For  the  calculation  of  sums  the  following  notations  will  be  used: 

REM     log(S-infinity  -  true  stress)  =  Y,  log(true  strain)  =  X, 

REM      IX  =  A#  ,  SY  =  B#,  IXY  =  C#,  IXA2  =  D#,  (IX)A2  =  A#A2  =  E#,  SYA2  =  F#, 

REM      (IY)A2  =  B#A2  =  G#  ,  R=  correlation 

REM 

CORR=  0  :  CORRMAX=  0  :  STP=  1     'CORR  and  CORRMAX  are  a  comparitive 

'     correlation  FOR  use  with  R  ;  STP 

'     is  a  step  size  which  will  decrease  as  we  narrow  in  on  S-infinity 

SINF*  SMAX +  1.11 

'     setting  S-infinity  equal  to  a  little  greater  than  SMAX  to 

'     prevent  the  possibility  of  logs  of  negative  numbers  and  we  know  S-infinity 

'     should  be  greater  than  SMAX 

PRINT  "Corr         S-infinity" 

QUE=0 
VOCESUM: 
A#=0  :  B#=  0  :  C#=  0  :  D#=  0  :  E#=  0  :  F#=  0  :  G#=  0  :  J=0  'initialize  variables  to  zero 

REM* 

REM**  Any  changes  to  this  area  must  also  be  changed  below  in  the  similar  areas 

REM******* 

FOR  l=  STARTPT  TO  STOPPT 

J=  J  +1        'counter 

A#=  A#+  STRAIN(I)        'sum  of  strains 
B0=  LOG(SINF  -  STRESS(I)) 

B#=  B#+  BO        'sum  of  the  dependent  variables 

C#=  C#+  (STRAIN(I)  *  BO)    'sum  of  the  products 

D#=  D#+  STRAIN(I)A2     'sum  of  the  squared  strain 

F#=  F#+  B0A2      'sum  of  the  squared  dependent  variable 
NEXT  I 

E#=  A#A2 

G#=  B#A2 
R=  ((J*C#)  -  (A#*B#))/SQR(((J*D#)-  E#)*((J*F#)-  G#))  'correlation 

itt******itir**  *********************  ************************ 

PRINT  USING  "+.#####     ###.##";  R.SINF 

IF  ABS(R)  >=  CORR  THEN     'have  found  a  correlation  which  is  higher 
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CORR=  ABS(R)  'save  that  correlation  as  the  new  one  to  beat 

IF  CORR  >=  CORRMAX  THEN 

CORRMAX=  CORR     'save  highest  correlation 

SINF0=  SINF  'save  corresponding  value  of  S-infinity 

END  IF 

SINF=  SINF+  STP        'increment  S-infinity  by  some  step  size 
GOTO  VOCESUM        '  return  to  summing  and  calc  routines 
ELSEIF  ABS(R)<CORR  THEN 
QUE=  QUE  +  1 
IFQUE<1  THEN 

SINF=SINF+STP 
GOTO  VOCESUM 
ELSEIF  QUE  >=  1  THEN 

IF  STP  <=  .1  THEN  GOTO  PUNCHOUT  'Escape  from  the  loops 

QUE=0 

SINF=  SINFO  -  STP  'go  back  two  full  steps  before  the  max  value  of  SINF 

STP=STP/1 0  'decrease  the  step  size 

CORR=0 
GOTO  VOCESUM 
END  IF 
END  IF 

i******************** ********** ********************* ******* 

PUNCHOUT: 

REM**  By  now  we  should  have  the  straightest  line  value  of  S-infinity  so  the 
REM**  slope  of  this  line  is  the  Voce  coefficient,  A,  we  will  call  it  M.  The 
REM**  intercept  is  So.  NEED  TO  RECALCULATE  THE  SUMS  OF  THE  LOG(...) 
REM**  BECAUSE  THE  MOST  RECENT  SET  IS  NOT  THE  STRAIGHT  LINE  SUM 
REM**  ANY  CHANGES  TO  THE  ABOVE  SECTION  MUST  ALSO  BE  REFLECTED  HERE 

A#=  0  :  B#=  0  :  C#=  0  :  D#=  0  :  E#=  0  :  F#=  0  :  G#=  0  :  J=0  :  SINF=  SINFO 

'reinitialize  variables 

FOR  l=  STARTPT  TO  STOPPT 

J=  J  +1         'counter 

A#=  A#+  STRAIN(I) 

B0=  LOG(SINF  -  STRESS(I)) 

B#=  B#+  BO 

C#=  C#+  (STRAIN(I)  *  BO) 

D#=  D#+  STRAIN(I)A2 

F#=  F#+  B0A2 
NEXT  I 

E#=  A#A2 

G#=  B#A2 
R=  ((J*C#)  -  (A#*B#))/SQR(((J*D#)-  E#)*((J*F#)-  G#))  'correlation 

REM     We  now  calculate  the  straight  line  slope  (M)  and  intercept  for  this  value  of 
REM     S-infinity.  The  intercept  is  the  quantity  (S-infinity  -  So)  which  we  will 
REM    call  SINTER,  SO  is  approximately  the  0.1  %  yield  strength. 
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M=  ((J*C#)-(A#*B#))/((J*D#)-E#) 

M1=  -1/M        This  is  required  because  Voce  takes  the  slope  of  the  line  where 
'  Strain  is  plotted  on  the  Y-axis  and  log(Sinf  -  S)  is  on  the  X-axis. ..the 
'  reciprocal  of  what  this  program  has  done  since  stress  is  usually  plotted 
■  on  the  Y-axis.  Where  this  'line'  crosses  the  stress  axis  is  the  location 
'  of  S-infinity.  This  value  of  the  (negative)  slope  is  called  the 
'  characteristic  strain. 

SINTER=  (B#/J)-(M*(A#/J)) 

S0=  -(EXP(SINTER)  -  SINFO) 

BEEP 
PRINT 
PRINT 

PRINT  "In  the  file  "  DAT$ 
PRINT  "Using  the  true  stress  and  true  strain:" 
PRINT 

PRINT  "The  Voce  equation:  S=  Sinf  -  (Sinf-So)e  A-(E/A)  ;  where  E  is  the  strain" 
PRINT  "Between  points  ";STARTPT;"  and  ";STOPPT;"  of  file  ";DAT$ 
PRINT  "Has  the  coefficients:       Si nf=  "SINFO"     So=  "SO"     A=  "M1 
PRINT  "With  a  correlation  of:  "R 

PRINT  "Necking  begins  at  point ";  MAXLODPT; "  where  the  load  is  ";  MAXLOD 
PRINT 
PRINT 
PRINT  "Do  you  want  a  printout  of  the  results  (Y  or  N):"; 

BEEP 

INPUT  ANS$ 

IF  ANS$=  "Y"  OR  ANS$=  "y"  THEN 
GOTO  PRINTOUT 

ELSE 
GOTO  RESTART 

END  IF 
PRINTOUT: 

LPRINT 

LPRINT 

LPRINT  "In  the  file  "  DAT$ 

LPRINT  "Using  the  true  stress  and  true  strain:" 

LPRINT 

LPRINT  "For  the  Voce  equation:  S=  Sinf  -  (Sinf-So)e  A-(E/A) ;  where  E  is  the  strain" 

LPRINT  "Between  points  ";STARTPT;"  and  ";STOPPT;"  of  file  ";DAT$ 

LPRINT  "Has  the  coefficients:         Sinf=  "SINFO"     So=  "SO"     A=  "M1 

LPRINT  "With  a  correlation  of:  "R 

LPRINT  "Necking  begins  at  point ";  MAXLODPT; "  where  the  load  is  ";  MAXLOD 

LPRINT 

LPRINT 

LPRINT  CHR$(1 2) 
RESTART: 

PRINT 

PRINT 

PRINT  "Do  you  want  another  run  of  this  data  set  (Y  or  N):" 
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PRINT  "Select  N  if  you  had  a  bad  data  input  and  Y  for  the  next  question"; 

BEEP 

INPUT  A$ 
IF  A$=  "Y"  OR  A$=  "y"  THEN  GOTO  BEGINNING 
PRINT 
PRINT  "Do  you  want  to  run  the  program  on  another  data  file  (Y  or  N):"; 

BEEP 

INPUT  A$ 
IF  A$=  "Y"  OR  A$=  "y"  THEN  GOTO  NEXTFILE 
END 
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APPENDIX  C 
BASIC  Computer  Program  For  Fitting  an  Exponential  Function 


REM**  THIS  PROGRAM  IS  TO  TAKE  TEMP/S0/S°°  DATA  FROM  A  DATA 

REM**  FILE,  PUT  IT  INTO  ARRAYS  AND  THEN  CALCULATE  THE  EQUATION 

REM**  COEFFICIENTS. 

REM** 

REM**  This  program  is  written  for  the  Apple  Macintosh  PC  in  Microsoft  BASIC 

REM**  Version  3.0 

i************************************************************ 

REM 

REM**  DIMENSIONING  THE  CALCULATION  ARRAY  AND  INPUT  OF  DATA 

DIM  TEMP(50) ,  S(50) 

NEXTFILE: 

PRINT 

PRINT 

PRINT  "Input  the  name  of  your  input  file  name  with  the  Temp(°C)/S°°  or  So" 

PRINT  "data  in  it.  The  data  should  be  in  the  same  order " 

PRINT 

PRINT  "INPUT  FILE  NAME:"; 

BEEP 

INPUT  DAT$ 
OPEN  DAT$  FOR  INPUT  AS  #1 
PRINT 
PRINT 
N=0 

REM**     LOOP  THROUGH  TO  STORE  TEMP/Soo  or  SO  DATA  IN  ARRAYS  FOR  FUTURE 
REM**     CALCULATIONS 

WHILE  NOT  EOF(1) 
N=N+1 

INPUT#1,TEMP(N),S(N) 
TEMP(N)=  TEMP(N)+273       'convert  °C  to  °K 

WEND 
CLOSE  #1 
REM 
AGAIN: 

PRINT  N;"  sets  of  data  have  been  read  from  file  ";  DAT$ 
PRINT 

PRINT  "Enter  the  range  of  points  you  desire  to  look  at  using  data  point  numbers" 
PRINT "  separated  by  a  comma."; 

BEEP 
INPUT  STARTPT.STOPPT 
PRINT 
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PRINT  "Sums  of  the  data  are  now  being  calculated  for  a  best  correlation  fit." 

PRINT 

PRINT 

REM     For  the  calculation  of  sums  the  following  notations  will  be  used: 

REM     stress  coefficients  =  Y,  Exp(Tc/T)  =  X, 

REM      IX  =  A#  ,  IY  =  B#,  IXY  =  C#,  XXA2  =  D#,  (IX)  A2  =  A#A2  =  E#,  IYA2  =  F#, 

REM      (IY)A2  =  B#A2  =  G#  ,  R=  correlation 

REM 

CORR=  0  :  CORRMAX=  0  :  STP=  1 0    'CORR  and  CORRMAX  are  a  co.mparitive 

'     correlation  FOR  use  with  R  ;  STP 

'PRINT  "Corr         TC" 

QUE=0 
VOCESUM: 
FOR  TC=  1  TO  300 

A#=0  :  B#=  0  :  C#=  0  :  D#=  0  :  E#=  0  :  F#=  0  :  G#=  0  :  J=0  'initialize  variables  to  zero 
REM******* 

REM**  Any  changes  to  this  area  must  also  be  changed  below  in  the  similar  areas 
REM******* 
FOR  l=  STARTPT  TO  STOPPT 

J=  J  +1         'counter 

A#=  A#+  EXP(TC/TEMP(I)) 

B#=  B#+  S(l) 

C#=  C#+  (S(l)  *  EXP(TC/TEMP(I))) 

D#=  D#+  (EXP(TC/TEMP(I)))A2 

F#=  F#+  (S(I))A2 
NEXT  I 

E#=  A#A2 

G#=  B#A2 
R=  ((J*C#)  -  (A#*B#))/SQR(ABS(((J*D#)-  E#)*((J*F#)-  G#)))  'correlation 

I-**********-***  A*************************  ***************** 

IF  ABS(R)  >=  CORR  THEN     'have  found  a  correlation  which  is  higher 
CORR=  ABS(R)  'save  that  correlation  as  the  new  one  to  beat 

IF  CORR  >=  CORRMAX  THEN 

CORRMAX=  CORR     'save  highest  correlation 
T0=  TC  'save  corresponding  value  of  TC 

END  IF 
END  IF 

NEXTTC 

i*********************#************ ************************ 

PUNCHOUT: 

REM**  By  now  we  should  have  the  straightest  line  value  of  TC  so  the 

REM**  slope  of  this  line  is  the  coefficient,  C1 ,  we  will  call  it  M.  The 

REM**  intercept  is  Co.  NEED  TO  RECALCULATE  THE  SUMS 

REM**  BECAUSE  THE  MOST  RECENT  SET  IS  NOT  THE  STRAIGHT  LINE  SUM 

REM**  ANY  CHANGES  TO  THE  ABOVE  SECTION  MUST  ALSO  BE  REFLECTED  HERE 

A#=  0:  B#=  0:  C#=  0:  D#=  0:  E#=  0:  F#=  0:  G#=  0:  J=0:  TC=  TO  'reinitialize  variables 
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FOR  1=  STARTPT  TO  STOPPT 

J=  J  +1         'counter 

A#=  A#+  EXP(TC/TEMP(I)) 

B#=  B#+  S(l) 

C#=  C#+  (S(l)  *  EXP(TC/TEMP(I))) 

D#=  D#+  (EXP(TC/TEMP(I)))A2 

F#=  F#+  (S(I))A2 
NEXT  I 

E#=  A#A2 

G#=  B#A2 
R=  ((J*C#)  -  (A#*B#))/SQR(((J*D#)-  E#)*((J*F#)-  G#))  'correlation 

REM     We  now  calculate  the  straight  line  slope  (M)  and  intercept  for  this 
REM    value  of  TC.  The  intercept  is  the  quantity  CO 

M=  ((J*C#)-(A#*B#))/((J*D#)-E#) 

C0=  (B#/J)-(M*(A#/J)) 

BEEP 
PRINT 
PRINT 

PRINT  "In  the  file  "  DAT$ 

PRINT  "For  an  equation  of  the  form:    S=  Co  +  [C1  *Exp(Tc/T)]" 
PRINT 

PRINT  "The  coefficients  are:" 
PRINT"  Tc="T0"     Co="C0"      C1="M 

PRINT  "With  a  correlation  of:  "R 
PRINT 
PRINT 
RESTART: 

PRINT 

PRINT 

PRINT  "Do  you  want  another  run  of  this  data  set  (Y  or  N):" 

PRINT  "Select  N  if  you  had  a  bad  data  input  and  Y  for  the  next  question"; 
BEEP 
INPUT  A$ 

IF  A$=  "Y"  OR  A$=  "y"  THEN  GOTO  AGAIN 

PRINT 

PRINT  "Do  you  want  to  run  the  program  on  another  data  file  (Y  or  N):"; 
BEEP 
INPUT  A$ 

IF  A$=  "Y"  OR  A$=  "y"  THEN  GOTO  NEXTFILE 
END 
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